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SUMMARY 


The basic goal of optimization of acquisition parameters is to have optimum 
resolution and signal to noise ratio in the acquired data. The aim and 
objective of the present book is to provide details about the optimization of 
on-land and offshore Seismic Data Acquisition Parameters. In this book I 
have tried my best to provide fundamentals behind the optimization so that 
user can quickly decide appropriate acquisition parameters needed for 
different situations. Accordingly the book has been divided into five 
chapters. 

Chapter-one gives an overview of the general technique of seismic data 
acquisition and common meaning of some of the terms being used in 
exploration seismic in on- land and offshore areas. 

Chapter-two deals with the optimization of Geometrical Parameters i.e. far- 
offset distance, near-offset distance, group interval, effective array length, 
foldage, line orientation and line length. To optimize the geometrical 
parameters, one must, first identify the main zone of interest and get the 
following information of the area 

• Maximum and minimum depths of the main reflector of 
interest in the area and corresponding average velocities. 

• Maximum depth of deepest reflector of interest, if any, and 
corresponding average velocity. 

• Minimum depth of shallowest reflector of interest, if any, and 
corresponding average velocity. 

• Minimum thickness of the bed to be resolved at main reflector 
level. 

• Maximum dip of the main reflector in the area. 

On the basis of above information, the geometrical parameters are decided 
in the following sequence. 

(i] Determination of shortest wavelength or highest frequency content 
of signal required to get proper resolution. 



[ii] Determination of far-offset distance, considering the requirement of 
velocity resolution, at the level of deepest reflector of interest, 
requirement of multiple attenuation and effect of NMO stretch at the 
main reflector level and to make appropriate choice. 

(hi] Determination of Near-offset distance, considering the minimum 
depth of shallowest reflector of interest and interference of ground 
roll with nearest channels. 

[iv] Determination of group-interval, considering the far offset distance 
and shortest wavelength component of signal and problem of spatial 
aliasing. 

[v] Determination of effective array length, considering the far offset 
distance, shortest wavelength component of signal and prominent 
noises wavelength in the area. 

[vi] Determination of line length, considering the migration aperture 
and foldage requirement at the end of the line. 

Chapter-three deals with the optimization of Source Parameters. For 
explosive source, shot hole depth, charge size and lithology of the medium 
surrounding the charge are three parameters which decide the source 
signature. After knowing the thickness of weathering layers in the area, 
records with different depths below the weathering layer with equal charge 
size is taken. Signal to noise ratio and peak frequency of the signal at the 
main zone of interest is computed. Amplitude spectra of different traces at 
main zone of interest are also compared for all the records. Record giving 
highest peak frequency and signal to noise ratio corresponds to optimum 
depth. Now keeping this depth fixed, charge size is varied and records are 
taken with different charges. Again signal to noise ratio, peak frequency and 
amplitude spectra are compared and the suitable one is selected on the basis 
of frequency content and signal to noise ratio. In the case of land air gun, 
number of air guns, number of pops and air gun patterns are the main 
parameters which decide signal to noise ratio. In the case of vibrator, 
sweep length end frequencies of the sweep and vibrator patterns are the 
main parameters to be optimized in the field. In offshore work, marine air 
gun is the widely used marine seismic source. Peak to peak amplitude, 
primary to bubble ratio [PBR] and bandwidth of the amplitude spectrum are 
the deciding parameters for proper selection of marine air gun source. 



Chapter-four deals with Noise Analysis and selection of appropriate 
Source/Receiver Array Parameters accordingly. Effective array length, 
element's weighting and spacing are the parameters to be optimized on the 
basis of desired seismic resolution as well as degree of attenuation of 
coherent noises. Effectiveness of the array is examined by fold-back 
experiment and f-k diagram. 

In Chapter-five, silent features of the Recording System & Quality Control 
have been discussed which are essential to be adhered in the field so that 
objective of the survey is achieved. 

Dr. R Prasad 
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CHAPTER-1 


SEISMIC DATA ACQUISITION: AN OVERVIEW 


1.0 INTRODUCTION 

T he principles of seismic data acquisition in on-land and offshore are 
similar. An artificial seismic wave (an elastic wave] is generated in the 
earth or in sea water by a suitable source. The elastic waves thus generated, 
travel into the earth and after reflections/refractions from different 
interfaces of sedimentary layers of earth crust the upcoming waves reach 
the ground surface. On the ground surface, in the case of land, several 
detectors called geophones are planted in a specified geometrical pattern, 
which convert the ground motion caused by reflected elastic waves into 
electrical voltages. In the case of sea, detectors are hydrophones placed 
below water surface, which convert pressure variation in water caused by 
reflected acoustic waves into electrical voltages. These electrical voltages 
are recorded in analog (true wave form] on photographic paper for quality 
control purpose and in digital form on a suitable device as a permanent 
record, for later use of processing and interpretation of these data. With the 
help of the time of arrival of reflected waves, its pulse shape, reflection 
pattern and its velocity in the earth medium, we can get the knowledge of 
the structure, stratigraphy and depth of the target horizon. A conceptual 
diagram for data acquisition is shown in Fig. 1.1 



Fig. 1.1: Conceptual diagram for data acquisition 
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1.1 SEISMIC WAVE 

When earth is excited with impulse of artificial seismic energy, several kinds 
of waves are generated. They may be divided into two categories (i] Body 
Wave and (ii) Surface Wave. Details for both kinds of waves have been 
discussed below. 

1.1.1 BODY WAVE 

This type of waves travels deep into the earth and therefore called body 
wave. There are two types of body waves. In the first type, the particle 
motion is along the direction of propagation and is called longitudinal or P- 
wave. In the second type, the particle motion is perpendicular to the 
direction of propagation and for this reason; this is called transverse wave 
or S wave. Velocity of P-wave is given by 

V 1= [(^+2p )/p] 1/2 

and velocity of S-wave is given by 

V 2 = (p/p ) 1/2 

Where p is density, X & p are Lame's constants, given by 
?i=aE / ((1 +g)(1-2g)) 

P = E/(2(1+g)) 

E is Young's modulus and g is Poisson's ratio 

1.1.2 SURFACE WAVE 

This type of waves are mainly confined to the upper most layer of the earth 
and hence called surface waves. There are two kinds of surface waves, the 
first is Rayleigh Wave and the second is Love Wave. In Rayleigh wave, 
particles oscillate in the vertical plane in an elliptic and retrograde 
movement. Amplitude of the wave decays exponentially as the depth 
increases. Vertical component of the motion is recorded by the vertical 
seismometers. Velocity (Vi ] of this wave may be given as 

Vi = 0.92(p/p) 1/2 = 0.92V 2 

V 2 being the velocity of transverse or S-Wave, 
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This wave shows dispersive character i.e. velocity varies with frequency. For 
high frequency i.e. for shorter wavelength, the wave does not penetrate 
deeply and obtains a velocity equal to 0.92 of the velocity of transverse wave 
in the upper bed; while for the surface wave of low frequency, and therefore 
greater penetration, the velocity will approach 0.92 of the velocity of the 
transverse wave in the deeper high velocity layer. 

In Love wave particle motion is perpendicular to the direction of 
propagation and in the horizontal plane. For this reason this is not recorded 
by vertical seismometer. 


1.2 SEISMIC NOISE 

In reflection seismology, we are interested in that P-wave which reaches to 
the surface as primary reflection, therefore called signal, and all other types 
of waves are unwanted and are called noises. The signal to noise ratio (S/N] 
is the ratio of the signal strength in a specified portion of the record to the 
noise strength in the same portion. Seismic noises can be broadly divided 
into two parts (i] Coherent noise [ii] Random noise 

1.2.1 COHERENT NOISE 

Coherent noise can be followed from trace to trace and have some common 
properties, at least up to few traces but incoherent noise is dissimilar on all 
traces and we cannot predict from one trace to other. The following are the 
common noises observed in seismic data. ( Fig.1.2] 

• Surface waves 

• Air waves 

• Short path multiples 

• Ghost 

• Near-Surface multiples 


Peg-leg multiples 
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• Long-path multiples 

(a] Simple multiples 

(b] Inter-bed multiples 

■ Diffraction 
• Refraction 

(a] Refracted P-wave 

[b] Refracted S-wave 

■ Reflected Refraction 

■ Mode converted waves (P-wave to S wave] 



SHORT PATH MULTIPLE LONG PATH MULTIPLE 

Fig. 1.2: Common Noises in Seismic Data 

1.2.2. NOISE OBSERVED IN MARINE DATA 


Following categories of noises are generally seen in marine seismic data. 
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• Direct arrivals of energy from the shot to the receiver through the 
water. 

• Multiple reflections or reverberations from the sea floor. 

• Refraction and multiple refraction from the sea floor. 

• Noise resulting from the mechanical motion of the cable in response 
to variable tugging force. The noises in this category are very low in 
frequency. 

• Noise scatterd from irregularities in the sea floor or irregularities in 
the very shallow earth beneath the sea floor. This category of noise 
some times are source of major problem in stack data at late record 
time. When there is sharp transition from water to a hard sea floor, 
there will be a lot of energy from the source trapped in the water 
layer. If sea floor is irregular, scattering of the trapped energy gives 
rise to strong events which exhibit moveout with velocity of sound in 
water. This moueout will vary from nearly linear to nearly hyperbolic. 
Scatterers that are not laterally displaced and are either forward of 
the vessel or beyond the tailbuoy will result in noise with linear 
moveout. Scatterers that are between the source and receivers and 
which are displaced a considerable distance, either laterally or 
vertically, will display a nearly hyperbolic moveout. 

• Swell noise is very common in rough sea conditions and is generated 
by wave action. This noise is dependent on cable depth, wave height 
and wavelength. Swell or wave noise increases with decreasing swell 
wavelength, it increases with increasing wave height and it decreases 
with increasing cable depth. The attenuaton and speed of the swell 
depend mainly on its wavelength. Long waves decay very slowly and 
can travel considerable distance from the area, where they 
originated. 


1.3 AIM OF SEISMIC DATA ACQUISITION 

The basic aim of data acquisition is to record signal of (i] good signal to 
Noise ratio (S/N] and (ii) broader bandwidth for better resolution to fulfil 
the geological objectives. These objectives can be achived by attenuation of 
noises and improvement of the signal strength as well as the frequency band 
width at the stage of data acquisition and data processing. The details have 
been discussed in the following sections. 
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1.3.1 SIGNAL TO NOISE RATIO 

Not all but some of the noises can be attenuated at the stage of data 
acquisition and some at the stage of data processing. Some part of the 
surface wave can be attenuated by frequency filtering; and some part with 
the help of proper source or/and receiver array pattern. Multiples can be 
attenuated at the stage of data processing if move-outs of primary and 
multiple reflections have enough separation. Strength of Signal can also be 
improved by selection of proper source parameters. 

The emphasis in recording should be to preserve high frequencies signal and 
suppress noises. The sampling rate and antialiasing filter should be 
adequate to record the desired frequencies. The length of receiver arrays 
should be small enough to prevent significant loss of high frequency signal 
due to intra-group move out and statics. However, the array length should 
not be too small, since small arrays are not as effective at suppressing 
random high frequency noise (wind noise] as large arrays. Finally, the 
source strength should be high enough to provide adequate signal level 
relative to noise within the desired frequency band. Unless the S/N of the 
field data is above some minimal value (say o.25] processing algorithms may 
have difficulty in resolving the signal. 

1.3.2 RESOLUTION 


(A) Vertical resolution 

Resolution refers to the minimum separation between two features such 
that we can tell that there are two separate features. If seismic wavelets 
were a spike, the resolution would not have been a problem. If the thickness 
of the thin layer is less than a certain minimum value, the top and bottom 
reflections from the thin layer cannot be seen separately. Raleigh criterion of 
resolution states that two events can be resolved if their separation is half 
cycle, since events are recorded in terms of two-way time, therefore real 
separation of the features must be quarter wave length. Thus resolvable 
limit is A/4 as seen in Fig.1.3 


Vertical Resolution = A/4 
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Fig. 1.3: Vertical Resolution 

Seismic wave velocities in the sub-surface range between 2000m/s and 
5000 m/s and generally increase with depth. On the other hand, the 
dominant frequency of seismic signal typically varies between 50 and 20 Hz 
and decreases with depth. Therefore typical seismic wavelength ranges from 
40 to 250 m and generally increases with depth. Therefore deeper feature 
must be thicker than the shallower features to be resolvable. In the case of 
fault, when the throw of the fault is equal to or greater than A/4 of dominant 
wavelength then the presence of fault can be inferred easily. 

(B) Lateral Resolution 

Lateral resolution refers to how close two reflecting points can be situated 
horizontally, so that they can be recognised as two separate points rather 
than one. Consider the spherical wavefront that strikes the horizontal plane 
reflector AA’ (Fig.1.4], Every point on this reflector can be considered as 
point diffracters. For a coincident source and receiver at the earth's surface 
'S', the energy from the sub-surface point 'O' arrives at 'S', if'Z' is the depth 
of the point diffractor 'O' and 'V' the velocity of the medium, then arrival 
time 

t = 2Z/V 

Now let the incident wavefront advance in depth by the amount A/4, energy 
from sub-surface location A, or A' will reach S at a time 

ti = 2(Z+ A /4)/V 

The energy from all the points within the circle of radius OA will arrive 
within time interval 
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ti — t = A/(2V) = T/2 (half wave period) 
and they interfere constructively. 

Two reflecting points that fall within this zone generally are considered 
indistiguishable as observed from the earth's surface. Spatial resolution is 
better understood in terms of diffractions. The reflecting disc AA' is called 
the first Fresnel Zone as shown in Fig. 1.4. The smaller the Freznel zone the 
easier is to differentiate between two reflection points Hence, the Freznel 
Zone width is a measure of lateral resolution. 



Fig. 1.4: Fresnel's Zone 

In the above figure, the radius of the first Freznel Zone is OA given by 
OA 2 = (Z+ A/4) 2 - Z 2 
=Z 2 (1+ A/4Z) 2 - Z 2 

= Z 2 (1+ A/2Z) - Z 2 since, Z >> A, (A/Z) 2 is neglected. 

= AZ/2 

OA = (Z A/2) 1 / 2 = {(1/2) (Vt/2) (V/f)} 1 / 2 =(1/2) V (t/f) 1 / 2 

Thus, the size of the zone depends on frequency, the higher the frequency 
the smaller the zone. Since, effective radius of the first Fresnel zone is half of 
the actual radius OA, we get 

Lateral Resolution = (1/4) V (t/f) 1 / 2 
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Where V = average velocity to the reflector, t= two way time and 
f= frequency. 


1.4 ACQUISITION TECHNIQUES 

Acquisition techniques may be described in three parts (i) Seismic Source, 
(ii) Seismic Receiver and (hi] Spread Geometry and Recording Instruments. 
Brief descriptions of these are given here and the details of optimization of 
related parameters will be discussed in the following chapters. 

1.4.1 SEISMIC SOURCE 

The first basic goal of data acquisition is to design a suitable seismic source 
(input function), so that it can bring more information about earth's interior. 
There are many types of Seismic Sources. Among them (a) Explosive (b) 
Vibrator and (c) Land Air Gun (d) Marine Air Gun are the four main sources 
of energy which are widely used. 

(A) Explosive 

In the case of explosive (dynamite), there are two factors which affect the 
input pulse shape very much. The first is proper depth of shot hole where 
the charge is placed and the second is right amount of charge size. Best 
information can be obtained if the input pulse has large amplitude and 
shortest possible time duration. Two types of explosives are used generally; 
(i) Gelatin dynamite (ii) Ammonium nitrate. The gelatin dynamite is the 
mixture of nitroglycerin and nitro cotton and an inert material which binds 
the mixture together and which can be used to vary the 'Strength' of 
explosive. The velocity of detonation is high for the explosives used in 
seismic work, around 6000 - 7000 m/sec; consequently, the Seismic pulses 
generated have very steep front in comparison with other energy sources. 
Electric blasting caps are used to initiate an explosion. These consist of 
small metal cylinders containing a resistance wire embedded in a powder 
charge. By means of two wires issuing from the end of the cap , a large 
current is passed through the resistance wire and the heat thus generated 
initiates the deflagration of the powder, which causes the explosion of an 
adjacent explosive in the cap. The current which is required for the blasting 
of cap, is derived from a blaster or encoder/decoder or syncrophones. These 
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are basically devices for charging a capacitor to a high voltage by means of 
either hand operated generator (blaster) or batteries (Encoder/Decoder & 
Syncrophones). Capacitor is discharged through the cap at the desired time. 
Biasing device also generates an electric pulse at the instant that the 
explosion begins. This time-break pulse fixes the instant of the explosion, 
(t=0). The time-break pulse is transmitted to the recording equipment by a 
pair of wires or radio where it is recorded along with the seismic data. 
Details about optimization of the source parameters will be discussed in 
Chapter-3. 

(B) Vibrators 

Vibrator is a surface energy source. Essentially all the surface energy 
sources have low energy compared to the explosive source. But vertical 
staking of several records can increase signal amplitude, all taken at one 
location without changing the receivers and by horizontal staking of several 
traces within a CDP gather after dynamic correction. Vibroseis sources pass 
oscillatory (not impulsive) energy of continuously varying frequency into 
the ground for several seconds. A control signal causes a vibrator to exert a 
variable pressure on a base plate pressed against the ground through 
hydraulic system. The pressure 'P' generally varies according to the relation, 

P = A(t) Sin 2 n t [fo + (df/dt) t ] 

df/dt is either positive (up - sweep) or negative (down- sweep) and is 
constant in case of ' linear sweep '. The amplitude A (t) is usually constant 
except during the initial and final 0.2 sec. or so when it increases from zero 
or decreases to zero. For linear sweep (df/dt=b), starting frequency is 'fo ' 
and ending frequency is (fo + bT), where T is the length of the sweep. Ideally 
the input to the ground is a copy of the pressure applied to the steel plate 
but in actual practice, crushing and compaction of the surface material 
causes the input to the ground varying non-linearly with the pressure 
exerted by the base-plate (pressure of the order of 200 Kg/Cm 2 ). This 
introduces harmonics, which are not present in the original input signal 
(sweep). Essentially the harmonics of the low frequency end of the sweep 
correlate with the proper fundamental frequency in the sweep signal and 
produce a ghost. 
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Fig. 1.5: Block Diagram of Vibroseis System 


Bandwidth, Central frequency, Sweep length and composites of the sweep 
decide the input vibroseis pulse. By cross-correlation of the input sweep and 
recorded seismic trace, the signal is compressed to a symmetrical pulse 
shape of large amplitude. Block diagram of vibroseis system is shown in Fig. 
1.5 

(C) Land Airgun 

Land air gun is another surface energy source. Although it is a low energy 
source, it is impulsive. Air compressed to very high pressure (2000 psi 
commonly] is discharged suddenly inside a bell shaped water chamber. 
Water chamber is sealed at the bottom by means of an elastomeric 
diaphragm and rest under its own weight in a circular base plate. The base 
plate is held in firm contact with the ground. Inside the bell, air gun is 
mounted. When a fire command is received, an electrically actuated valve 
triggers the air gun. The explosive release of high-pressure air into the water 
causes instantaneous expansion of the diaphragm creating a pressure wave 
into the earth and thus the sharply defined source signature is produced. 
The explosive release of the air occurs in 1-4 ms while the entire discharge 
cycle requires 25-40 ms, these are common values. Block diagram of land 
airgun system is shown in Fig. 1.6. Signal to Noise ratio depends on the 
number of air guns used at a time, number of pops at a place and pattern of 
airgun used. 
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Fig. 1.6: Block Diagram of Land Air Gun system 


(D) Marine Air gun 

The basic function of an under water source is to introduce a sudden 
pressure impulse into water. The most popular seismic source for marine 
work utilizes groups of guns. The guns are steel chambers, which are filled 
with high-pressure air that is released into water to generate the seismic 
signal. The air exhausted into water generates a spherical air bubble. Air 
bubble starts expanding and contracting due to elastic property of water. 
This oscillating bubble also creates seismic waves on each oscillation. After 
few oscillations the bubble reaches the water surface and releases the air. 
Period of bubble oscillation T depends on chamber volume of the gun as 
given below. 

T= 1.14 p 1/2 P 0 ' 5/6 Q 1/3 

T= Period of bubble oscillation 

p = density of water 

Po= Surrounding pressure 

Q= Potential energy inside the bubble, which is related to pressure 
and volume of air. 

Therefore, the bigger the chamber the longer is the period of oscillation. A 
typical gun signature is shown in Fig.1.7 
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Thus the seismic signal from an airgun includes the initial pulse plus the tail 
of subsequent low frequency pulses produced due to bubble oscillations. To 
minimize the effect of pulse repetition, several guns of different volumes are 
fired in synchronization. All initial pulses add in phase but bubbles having 
different periods of oscillation do not add in phase. This way the influence of 
the bubble oscillations is considerably reduced. A tuned airgun array is one 
that has been designed such that the undesirable bubble oscillations are 
minimized and dominant seismic signal is from the initial pulse when the 
high-pressure air is injected into the water. Peak to peak amplitude, primary 
to bubble ratio (PBR] and bandwidth of signal are main parameters to be 
decided. All these factors will be discussed in chapter-3. 



1.4.1 SEISMIC RECEIVER 

Detector of seismic wave is a simple device for the conversion of mechanical 
energy into electrical energy. This works on the principle of electromagnetic 
induction or peizo-electricity. 

(A) Geophone 

The detector used on land is called geophone, witch consists of a permanent 
magnet and a coil. A permanent magnet is freely suspended by a suspension 
spring with the casing and coil is firmly attached with the casing. On arrival 
of the wave, the earth's particles begin to vibrate and so the coil. Therefore 
there is a change of magnetic flux linked with the coil and hence e.m.f. is 
induced in the coil, which is sent to recording instrument through cable. This 
type of detectors is velocity sensitive. Our aim is to transform all motions of 
the earth into proportional electrical signal i.e. to make the geophone 
respond approximately equal for all frequencies above the natural 
frequency. 
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Performance characteristics of a geophone include (a] Natural Frequency 
(b] Sensitivity (c] Damping factor and (d] Harmonic distortion 

Natural Frequency 

Natural frequency Tn' is the frequency at which the undamped geophone 
oscillates when freely excited, natural frequency is a function of the spring 
constant K and the mass M. 

f N = (1/2 n) (K/M)i/2 

When geophone is excited at this frequency by earth motion resonance 
occurs. 

Sensitivity 

Sensitivity 'G' is closely related to magnetic field strength 'B', number of 
turns in the coil 'N' and the Cosine of the angle ' 0' between coil axis & the 
direction of magnetic field. 

G = BNC CosO 

Where, C is the geometrical constant of the geophone. Typical value for 
geophone sensitivity is 0.135 to 0.150 Volt/cm/sec. 

Damping Factor 

For damped system, the amplitude response at resonant frequency depends 
upon the damping factor. Damping flattens the amplitude response and 
reduces phase distortion. However it also reduces the sensitivity. Since 
oscillatory nature of geophone is undesirable, commonly 70% of critical 
damping factor is preferable. Damping factor ' [3' is related to the successive 
amplitude of damped oscillation (Fig. 1.8] 



P = Ln (A 1 /A 2 )/ [n 2 + {Ln (A 1 /A 2 )} 2 ] 
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Where, Ai & A 2 are amplitudes of consecutive half cycles. 

Harmonic Distortion 

Harmonic distortion results from the non-linear characteristics of the 
geophone or divergence of stress-strain relationship from Hook's Law. A 
typical tolerance value for total harmonic distortion is 0.2% or less at 
resonant frequency. 

It is desirable to use geophone in the frequency range above the natural 
frequency because amplitude response is flat and has little phase distortion. 
Therefore, design trends were to build geophones of lower natural 
frequency. However current attitudes favor using high frequency geophones 
because geophone response can be used as a low cut filter and because the 
operating characteristic of a high frequency geophone with high frequency 
data are better than those of a low frequency geophones recording the same 
ground motion. When used as low cut filter geophone reduces the spectral 
misbalancing caused by low pass earth filter and makes the recording of 
high frequency data possible by attenuating the low frequency component of 
signal. So that it cannot saturate the gain of recording system. 

The geophone ground coupled system can be represented as a damped 
harmonic oscillator. The characteristics of this coupling system depend upon 
the weight of the geophone, the effective area of contact, amplitude of 
ground motion and the elastic moduli of the ground. Coupling differences 
among geophones cause distortions in the recorded wavelet from trace to 
trace. These difference can be due to variation in plantation or change in the 
elastic parameters of the ground. One way to increase the effective area of 
contact is to increase the size of the spike or base plate of geophone. 
Geophone output decreases as the angle of coil axis with the vertical 
increases. 

(B) Hydrophone 

The detector used in sea water is called hydrophone which is piezoelectric 
crystal and is pressure sensitive. Certain crystalline materials develop an 
electric potential when they undergo a physical deformation. These 
materials are said to be piezoelectric (quartz, tourmaline, etc.]. Hydrophone 
is lowered in seawater. On arrival of the wave, the pressure on the crystal 
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surface varies and therefore electrical voltage is developed across the 
surface, which is sent to instrument through cables for recording on the 
magnetic tape or any suitable recording media. Hydrophones are equivalent 
to a pressure-controlled source in series with a capacitor. Following are the 
characteristics of Benthos Geopint hydrophone: 

Voltage sensitivity = 20V/Bar ± 3.5V/Bar 

Charge sensitivity = 320 nC/Bar ± 40 nC/Bar 

Capacitance = 16 nF ± 10 % 

For array of six hydrophones in parallel (6.25 m group interval] 

Charge sensitivity = 1920 nC/Bar 
Capacitance = 96 nF 

For array of 12 hydrophones in parallel (12.5 m group interval] 

Charge sensitivity = 3840 nC/Bar 
Capacitance = 192 nF 

1.4.3 SPREAD GEOMETRY AND RECORDING INSTRUMENT 

To acquire seismic data over an area, several lines (or profiles] are run 
through the area of investigation in a grid pattern. Direction of the lines 
depends on the dip and strike direction of the interface to be mapped and 
length & spacing of the lines depend on the size of the structures to be 
delineated. Along each line at regular interval, seismic pulses are 
transmitted into the earth and the detectors, also located at regular interval 
along the line, receive reflected pulses. All the detectors are connected to 
recording instrument through a multi-wire cable or radio wireless. These 
pulses are converted into electrical pulses by the detectors and recorded on 
magnetic tape or on the paper or any recording media by the recording unit. 
Distance between two consecutive source/shot points is called 'Source/Shot 
Point Interval' and the distance between the two consecutive receiver points 
is called "Group Interval”. The horizontal distance between source point and 
receiver point is called "Offset distance " (X] and travel time of the pulse 
from source to receiver via reflector is called 'two-way reflection time' (t]. 
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(A) 2D-Land-Seismic Data Acquisition 

In Figure 1 . 9 , there are four active channels/receivers connected to the 
recording unit at a time. Energy from the source Si is received 
simultaneously by the receivers Ri, R2, R3, R4 after reflection from the mid 
points Pi, P2, P3 and P4 between source & receiver of a plane horizontal 
reflector. Now source point is moved from Si to S2 and receivers from R2 to 
Rs are connected with the instrument, which get energy reflected from 
points P3, P4, Ps and ?(,. Similarly, the receivers R3 to R 6 receive energy from 
the source S3 after reflection from the points P5 to Pa This process continues 
till the end of the line is reached. In this case we see that each reflection 
point (except Pi and P2J is mapped twice through different length of ray 
paths. Thus we get 'two fold data' and this method of data acquisition is 
called C.D.P. (Common Depth Point] method, because the reflection points 
are common for several source and receiver positions. Common depth 
points and common mid points coincide in the case of plane horizontal 
reflector but for inclined reflector may have different locations. CMP interval 
is half of the receiver group interval. Fold of the data is the number of times 
the reflected pulse (event] are recorded from the same reflection point 
through different combination of source and receiver having different 
lengths of ray paths. Thus 

Fold = Zi (Number of receiver groups X group interval)/ shot point 
interval 


Above method holds for any number of channels. 



Fig. 1 . 9 : Multifold 2 D data acquisition in land 
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As the number of active receivers is increased, foldage of data is also 
increased. In multi-channel recording, distance of nearest receiver from the 
source point is called 'near trace offset' and distance of farthest receiver is 
called 'far trace offset'. Length of time up to which events (data] are 
recorded on the tape or paper is called 'record length'. 

End-On and Split- Spread CDP Technique 

When all the receivers are located on one side of the source point, it is called 
'end-on' technique as shown above in Fig. 1.9. When receivers are located 
on either side of the source point, it is called ' split spread ' technique. If 
number of receivers is equal on either side of source point, it is called 
'symmetric split spread ' and if unequal it is called 'asymmetric split spread'. 
The goal of CDP technique is: 

• To attenuate multiple reflection by out of phase addition in data 
processing. 

• To attenuate random noise by out of phase addition in processing. 

• To provide statistical data for velocity analysis. 

• To provide statistical data for static correction. 

The split spread is the best when data enhancement by the cancellation of 
random noise is primary objective. The split spread does not have to be 
symmetric about the shot point. The number of channels on one side should 
be established on the basis of the desired far-offset distance. The remaining 
number of channels can then be placed on the opposite side of the shot 
point. Split spread shooting provides a forward and reverse profile. This 
may be useful in data processing with dipping reflectors. End on shooting 
has the advantage of better multiple cancellation in CDP stack, because all 
channels have independent ray path geometry. The ray paths in split spread 
are not all independent because the traces in each CDP gather can be paired, 
each pair differing only in having shot and detector interchanged. Thus the 
resultant doubling of multiplicity should not improve attenuation of shot 
generated noise such as multiple reflections and ground rolls. A weighted 
stack may be desirable in processing of split-spread data to achieve 
satisfactory multiple attenuation. The asymmetric split spread is particularly 
useful when both shallow and deep horizons are targets. It provides 
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increased shallow multiplicity compared with end on shooting and longer 
far-offset channels than a symmetrical split spread. 

The data acquired as discussed above is called 2-D data, because we assume 
that the reflection points lie in the vertical plane (x, t] passing through the 
profile; direction of the profile is X-axis and two-way reflection time/vertical 
direction is t-axis. In this case, data acquired along the profile give the idea 
of sub-surface geology along the line of intersection of reflector with vertical 
plane passing through the profile. We cannot say anything about other 
directions. In the area of complex geology the ray reflected from a point, 
which is not in the vertical plane may also reach the receiver located in the 
vertical plane and the ray reflected from the point in the vertical plane may 
not reach the receivers located in that plane. This will cause distortion of 
subsurface picture. To overcome this effect in the area of complex geology 
3-D data acquisition is carried out. 

(B) 3D- Land Seismic Data 

All the principle of 2D land data acquisition are also applicable to 
acquisition of 3D land data. Land 3D acquisition commonly is carried out by 
swath shooting in which receiver cables are laid out in parallel lines (in line 
direction] and shots are positioned in a perpendicular direction (cross-line 
direction]. In Fig. 1.10, geometry of swath shooting with four receiver lines 
is shown. In the above case, as the first shot is fired, data not only from one 
CMP line is received but four CMP lines over an area abed of the reflector is 
received in terms of two way time. Thus we are mapping in (x, y, t] plane, 
where X-axis is the direction of receiver line, Y-axis is the direction 
perpendicular to the receiver line and t-axis is the depth in time. 

As the first shot line is completed, the receiver cables are moved up along 
the swath a number of stations equivalent to the shot-line spacing and 
shooting is repeated. Once one swath is completed, another swath parallel to 
it is recorded; this procedure is repeated over the entire survey area. 
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Binning 

In order to see the areal distribution of foldage in sub-surface, the survey 
area is sub-divided in rectangular grids. These grids are called bins. Binning 
programe counts the number of CMP in each bin. The process of putting 
traces into bins is called binning. Binning of the data is done in a cell or bin 
of following size: 

Minimum inline bin size = half the group interval 
Cross-line bin size = sub-surface CMP line interval 

CMP's falling in one bin are stacked together in processing and placed in the 
centre of the bin. 

Telemetry Recording System 

In recent land survey, large numbers of channels are used particularly for 
3D data acquisition. If numbers of channels are very large, it is not possible 
to acquire data using conventional cables in which each channel has pair of 
conducting wires for data transmission. As the number of channels 
increases cable become more & more bulky and total length of the cable 
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increases to such an extent that equivalent damping of the geophone is 
appreciably changed at far off distances, which degrades the data. To 
overcome this difficulty Telemetry data acquisition is carried out whenever 
large number of channels is used. In this type of data acquisition, a group of 
channels (generally 8 or 12] are connected through conventional cable to a 
unit called Data Acquisition Unit (DAU) or Array Terminal Unit (ATU].This 
name is different for different manufacturer. DAUs or ATUs are kept on the 
line at equal interval. Each DAU or ATU is connected to Central Recording 
unit or central control unit, through a separate cable or antenna. Sampled 
data is directly sent to central control unit serially by each DAU or ATU for 
formatting and recording on tape. Block diagram of data telemetry 
acquisition system is shown in Fig. 1.11 
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Fig. 1.11: Block diagram of Telemetry System 

(C) 2D-Marine Seismic Data 

In off-shore, the receiver cable (called streamer] and source arrays along 
with other in sea equipment attached with these are towed by the ship. First, 
steamer is lowered into water and positioned in the center of deck. Then gun 
arrays are lowered. One sub-surface CMP line is generated just below the 
sail line after complition of firing on the particular sail line. 2D Marine 
geometry is shown in Fig.1.12. Main elements of in-sea and on-board 
equipments, timing sequence and details of other considerations are 
discussed below. 
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Gun Array 

Generally two gun arrays one on either side of streamer or four gun arrays 
two on either side of the streamer depending on the requirement of total 
gun volumes are lowered. All the gun arrays fire simultaneously within a 
predefined time window and constitute a single source. Center of the gun 
arrays determines the position of the source. 

Active Section 

Active sections are the sections consisting of the receivers. Receiver consists 
of an array of hydrophones. Hydrophone measures the pressure variations 
with the aid of piezoelectric material, which generates a voltage upon 
deformation. Geometry of hydrophone array has effects on the frequency 
response. A streamer in water will be subject to movement in all directions 
due to currents and wave action. Therefore, two piezoelectric elements in 
one hydrophone are connected and polarized so that voltages due to 
pressure waves add and voltage due to directional acceleration will cancel. 
The hydrophones and other elements within the streamer are surrounded 
by oil, which maintains propper buoyancy of the cable in water. 



Fig. 1.12: 2D Marine data acquisition 

Lead-In-Cable 


The lead- in cable is the first front end part of the in water elements of a 
streamer. It runs from the reel on the vessel to the front end of the streamer. 
The lead- in is composed of a hard plastic skin covered with armoured 
shield, which protects the carrier wire inside. Because of its stiffness and 
armour, the lead- in will not easily be damaged by turbulance from the ship’s 
motion or by the pressure waves from the gun arrays. Towards ship end 
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lead- in is connected to intrument through another cable, called deck cable. 

Stretch-Section 

Between the lead in and active streamer sections a few strech sections are 
connected. Strech sections contain nylone ropes instead of steel ropes as in 
live sections. Nylone ropes make it more elastic. The function of strech 
section is to absorb the jerks caused by the ship. Another strech section is 
also connected at the far end of the streamer to absorb the jerks caused by 
tailbuoy. 

Electronic Module or Bubble 

In all recent marine data acquisition system using digital streamer, 
electronic modules also called bubbles are connected between sections of 
cable, which interface the hydrophones to the onboard acquisition system. 
Bubbles, upon receiving command from the onboard acquisition system, 
start sampling and digitizing the signals from all the hydrophone groups 
almost simultaneously. Digitized signals are then transmitted via 
transmission lines in the cable on a time-multiplexed basis by all bubbles to 
the onboard recording system. 

Depth-Contoller Bird. Depth Sensor and Compass Bird 

At regular interval through out the cable several depth controller birds are 
connected. These birds are programmed to keep the streamer at a preset 
depth, which they accomplish by using their wing angles. They are 
connected with the help of coders which can be attached with the streamer 
at desired location. Birds are inductively coupled with a coil in side the 
sectin through which they communicate with the onboard depth controller 
unit. Birds has also got a depth sensor, which measure the depth.The same 
depth cintroller bird may also have heading sensor or there may be separate 
birds having heading sensor, and they are attached to the cable at regular 
interval in the same way as depth contoller birds are. They measure 
direction of the cable at sensor position from magnetic north. Therefore 
correction for magnetic declination is to be incorporated to get direction 
with respect to true north. The compass data is input of a streamer shape 
algorithm by whivh actual shape of the sreamer is detemined using least 
square polynomial curve fitting method. Since the earth magnetic field is 
subject to various changes, the accurate correction for magnetic declination 
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is difficult. For this reason, the accoustic location systems are becoming 
more popular since they are able to determine the shape of the streamer at 
some locations through which the compass can be caliberated. 

Tailbuoy 

The tailbuoy is the end component of stremer having GPS unit for sattelite 
positiong. Apart from its role in position determination of streamer it also 
marks the end of the streamer. Since the tailbuoy can easily be located from 
the vessel, this is especially useful if the streamer breaks. 

Fathometer 

A fathometer is a device used for measuring the water depth with high 
accuracy. A tranducer converts electrical voltage into mechanical vibrations 
which are sent through water, these high frequency vibrations got reflected 
from the sea bottom and the echoes are received by tranducer, which then 
acts as a microphone. The received sound waves are converted to varying 
voltages which is ampified and transferred to the recording chart. Using the 
time of arrival of reflected signal and assumed velocity of sound in water, 
depth of sea bed is calculated. The collected data is digitized and sent to the 
navigation computer for recording on tape . 

Acquisition System 

This system basically consists of two units. First is seismic intrface unit. This 
unit is the shipborne interface between the digital streamers and onboard 
recording system. After the system parameters have been entered and the 
streamer has been initialised, the interface unit is ready for collecting data 
from in sea electronic modules or bubbles. This unit also performs various 
tests i.e. hydrophone leakage test, amplifire impulse response test, RMS 
noise test and DC offset test on electronic modules. The second sub system is 
the writing system, which writes data on cartridges in SEG-D format. Control 
over the two subsystems is provided by the system controller, which is 
computer based unit and has a high speed serial link to both the subsystems. 
A mouse driven user interface to the system controller handles system 
operation and parameter entry. Real time cable data and other QC displays 
can also be viewed. Acquisition system is also connected to the positioning 
system. 
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Source Control System 

This system performs all control, timing and firing functions of the source 
array. This system performs statitical analysis and shows graphic display of 
the gun performances, like no fire, autofire and fire out of specified window. 
The airgun timing controller is based on the principle of electromagnetic 
induction. A coil is located in the head of each gun and the top portion of the 
moving shuttle has several small magnets embedded into it. When the gun is 
fired, the magnets enter the coil thus inducing an electrical current which is 
recognised by the source controller. This way the controller determines the 
actual firing times of the guns with a resolution of 0.1 ms. With each shot 
point a graph is generated displaying gun firing time for every air gun. 

Bird Control System 

Bird Controller module allows control over the birds. Depth commands can 
be given to individual bird through this module and also battery voltages can 
be checked. Since birds contain a depth and heading sensor, the system can 
generate a view of the depth and shape of the streamer. 

Gyrocompass 

The gyrocompass consists of a balanced spinning wheel of which spin axis 
aligns with the earth’s axis of rotation. This is used to determine the 
direction of the true north. Dynamic errors can occur due to the roll of the 
vessel which can be minimised by mounting the gyro low on the vessel and 
at the centre of the roll. Further more gyro is placed in alignement with 
respect to the fore and aft line of the vessel. 

Positioning System 

Positioning system consists of two subsystem. First is absolute positioning 
sustem, which makes use of GPS and DGPS to determine the accurate real 
time position of the vessel's navigation antanna also called Integrated 
Navigation System (INS], Second is relative positioning system which makes 
use of the accoustic positionig system and the compass, also called Seismic 
Integrated Positioing System (SIPS] which determines the real time position 
of receivers and sources with respect to vessel's navigation antenna. SIPS is 
mainly used in 3D-marine data acquisition in the case of multi streamers and 
multi sources. 
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Global Positioning System (GPS] is based on technique of code and carrier 
tracking of sattelite signal transmitted by various sattelites. There are 
several sattelites flying in different orbits. They make one revolution in 12 
hours at an altidude of 20200 km. The information which they transmit 
consists details of the orbit they fly in and this enables the user to determine 
its position relative to the sattelite. In order to make sure that sattelites 
transmit the correct orbital information, ground stations measure the 
trajectores in which the satellites fly. These measurements are used to 
calculate the orbit constants after which this information is stored into 
sattelite by radio. In one of the method, travel time of the sattelite signal is 
measured and multiplying this by the speed of light the range of each 
sattelite is calculated by onboard receiver. In order to measure the travel 
time, the receiver generates signal identical to the transmitted signal and by 
cross-correlation of these two signals, the time delay of the source signal can 
be determined. The resulting range is called pseudo range and accuracy is 
about 300 m. There are four unknowns i.e. 3 position coordinares [x,y,z] of 
navigation antenna and fourth is receiver clock offset from the sattelite 
clock. Therefore, atleast 4 satellites data are needed at any time to solve 4 
range equations. At the same time, reference stations near to survey area 
also receive the same data from the sattelite. The reference stations then 
calculate the positions relative to the satellites. Since the exact position of 
the reference station is known, the corrections needed for accurate sattelite 
positiuoning can be determined. This is called differential correction. Once 
the corrections are calculated they are sent to the user by radio. The user 
applies these corrections to his own calculation, which results in 4 accurate 
ranges from which the position can be determined with an accuracy of 5 m. 
various QC calculations are also performed by the system to evaluate the 
quality of positioning data. 

Timing Sequence 

Shot point distance along the line is calculated by INS system, which also 
triggers the other system when shot point is encountered. It does so by 
sending a closure to the sequence timer. This sends a closure to the 
acquisition system as a system start. The acquisition system is then armed 
and prepared to receive a field time break [FTB] at which point the 
recording starts. The sequence timer sends a timing closure to the source 
controller, say 100 ms after it received the closure from INS. Less than say 
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20 ms later, the source controller sends a FTB to the acquisition system and 
data acquisition begins. The gun array is fired say 51 ms after the fire 
closure have been received by the source controller. 

Streamer Feathering 



The receiver cable is subject to a certain amount of sideways drift from the 
ideal sail line direction, which is called feathering. The angle between the 
actual cable position and the sail line direction is called feathering angle. 
Assuming a simple cable shape as shown in Fig.1.13, 10°feathering angle 
will cause around 265m shift in mid-point and 530m at tail end of receiver 
for 3 km cable length. 

Shot and Receiver Ghost 

Let us consider an airgun at a depth of d g below the water surface. When the 
gun is fired, seismic energy propagates outward in all directions. The 
reflection coefficient of water-air interface is -1 approximately. Therefore 
the ray that travels vertically upward is get reflected from the water-air 
interface with polarity reversal. It then travels vertically downward. This 
wave is called source ghost. Path difference between the source ghost and 
direct wave is given by (Fig. 1.14] 

A =2d g 

and total phase difference 

O g = II + 4ndgA < n was added due to phase reversal of ghost. 

We can represent direct and the ghost wave by two vectors separated by the 
angle <|)g as shown in Fig. 1.16. If amplitude of the individual wave is A, then 
the magnitude of the sum of the two vectors would be 

Ai = | 2A Cos (Og/2] | 

Ai=| 2A Cos (n/2 + (2n d g A]| 
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= | 2A Sin (211 d g /A)| 

Resultant amplitude is zero when 

2 d g /A = n = 0, 1 , 2, 3,. 

So, there will be first notch at 
2 d g /A =1 or A =2 d g or f = V/ (2 d g ) 
Where, V is velocity in water. 



Fig. 1.14: Fig. 1.15: Fig. 1.16: 

Source Ghost Receiver Ghost Vector Sum 

There will be first maximum, when 
2 d g /A =1/2 or A = 4d g or f= V/ (4d g ) 

For d g = 5m and V= 1500 m/s 

Frequency of first notch = 150 Hz. 

Frequency of first maximum = 75 Hz. 

Similarly at the receiver, there is direct arrival of reflections from reflector 
below and another arrival of reflections from the water-air interface above 
which is called receiver ghost as shown in Fig.1.15. 

The equations in this case are identical to those for the source ghost except 
the cable depth d c and phase <()c replacing the gun depth d g and phase i|) g . 

The resultant amplitude will be 

A 2 = | 2A Sin (211 d c /A]| 
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The first notch will occur at 
X = 2 d c or f = V/ (2d c ] 

For d c = 10 m, V= 1500m/s 
Frequency of the first notch =75 Hz. 

Frequency of the first maximum = 37 Hz. 

The effect of the above two ghosts are multiplicative, therefore for a smooth 
sea the total response would be 

A = Ai. A 2 

A = 4A 2 Sin (211 dJX) Sin (211 d c /k) 
with a phase angle <)) = (t))g/2 +<(>c/2] 

If guns and cables are operated at relatevely shallow depth, the response 
would be improved at relatively higher frequencies but would be reduced at 
lower frequencies. Conversely, if the guns and cables are operated at greater 
depths the response at lower frequencies would be improved but reduced at 
higher frequencies. To get wider bandwidth streamer depth should be kept 
different from the source depth, so that the corresponding notches occur at 
different frequencies. 

Streamer balancing 

Currents push the streamer out of the line direction if they come from the 
beam and creating feathering. On the other hand currents parallel to the line 
affect the stremer balancing. If the vessel is moving up the current the 
relative speed is high. Therefore the birds have enough flow around their 
wings to keep the cable at the right depth. Only the front end of the streamer 
may rise a little because of higher resistance of the cable. When the vessel is 
moving down the current the relative speed of the streamer with respect to 
water decreases, which results in less effective depth control by the birds. 
Furthermore the streamer is slack and thus some parts of the streamer may 
sink while other parts rise. The front end of streamer generally sinks 
because of greater weight of lead-in cable. 
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(D) 3D-MARINE DATA ACQUISITION 

A marine 3D survey involves shooting a number of closely spaced parallel 
2D lines, this is known as line shooting. Therefore, all 2D data acquisition 
equipment are used as discussed above in the case of 3D also. In addition 
one other equipment is used, which is called binning and display system. If 
multi-streamer and multi-source are used then apart from the above, 
equipment called seismic Integrated Positioning System (SIPS] is also used. 
For a particular case of 3 streamers and 2 sources (4 gun arrays], first of all 
outer two streamers are lowered and are dragged outside with the aid of 
paravanes. In this way desired separation between the streamers is created 
and maintained. Then the third streamer is lowered and positioned in the 
center. When the streamer cables are in position, 4 gun arrays are lowered 
and positioned in pair on either side of the central streamer. Pair of arrays 
on right side (starboard side] of central streamer fire first and constitute 
source-1. Similarly the pair of arrays on the left side (port side], which fire 
later constitute source-2. Centers of the respective pairs determine the 
position of two sources. When source-1 makes fire CMP lines 1, 3 and 5 are 
generated. CMP lines 2, 4 and 6 are generated when souce-2 makes fire. If 
streamer separation is 200 m and source separation is 100 m, then CMP line 
separation will be 50 m. 

Real Time Binning 

In order to show the areal distribution of the foldage in sub-surface, the 
survey area is subdivided in rectangular grids. The rectangular grids are 
called bins or cells. Binning programme counts the number of CMP in each 
bin and result is shown as a color-coded coverage map. The process of 
putting traces into bins is called binning. 

To determine the content of each bin the following basics rules are used: 

• Only one reflection point will be accepted for each offset. 

• If no prime line trace is available, a trace from adjacent line may be 
accepted. 

• If more traces from adjacent lines are available, the one closest to the 
center of the bin is accepted. 

CMP's falling in one bin are stacked together and placed in the center of the 
bin if feathering is small. The large feathering spreads the reflection points 
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all over the bin and thus a weighted mean of these locations has to be 
calculated. The resulting center of gravity is accepted as the reflection point 
of the summed trace. If two adjacent lines both suffer from large feather, the 
missing traces can be compensated by certain percentage of bin overlap in 
cross-line direction. Both bins share traces falling within the overlap zone. 
Each bin should contain at least a certain minimum number of traces. For 
binning purposes the whole streamer is generally divided into four equal 
segments. Each segment is assigned a minimum fold. They are generally 
90% for the near segment, 80% for the near mid segment, 75% for the far 
mid segment and 60% for the farthest segment. Overall, one bin should thus 
have minimum of 76% of total fold. 

Single Streamer Single Source 

If only one streamer and one source are used, one CMP line is generated for 
each pass of the vessel as shown in Fig.1.17. In this case 

Cross line bin size = sail line spacing 

Minimum inline bin size= % receiver group interval 

Fold = % (number of receiver groups X group interval] /shot interval 

Bin fold = (number of receiver groups X inline bin size] / shot interval 
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Fig. 1.17: Single Streamer Single Source (without feathering] 

Multi Streamer Multi Source 

If multi-streamers and multi-sources are used, then it is possible to acquire 
several CMP lines in one pass of the vessel. In Fig. 1.18 Geometry of 3- 
streamers with dual sources is shown and 6 CMP lines are recorded during 
each pass of the vessel. In this case 
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Cross line bin size = % source separation 
Minimum inline bin size= Vi receiver group interval 


Fold = % (Number of receiver groups X Group interval] 

(Shot interval X Number of sources] 

Bin fold = (Number of receiver groups X Inline bin size] 

(Shot interval X Number of sources] 

Effect of Feathering on CMP Distributions 

The receiver cable is subject to a certain amount of sideways drift from the 
ideal sail line direction, which is called feathering. The angle between the 
actual cable position and the sail line direction is called feathering angle. The 
maximum feathering angle for ideal coverage can be determined by the 
angle 

O max = Sin - 1 (2.16AX l /L) 

AXl = Line spacing 


L = Streamer length 



Fig. 1.18:3 Streamers 2 Sources with feathering 

If feather angle is higher than this value, it ih not necessarily a problem if 
adjacent lines have also similar feather angles. 

Because of significant variations in cable shape due to varying feathering 
angle during the course of recording, the mid point distribution in the cross- 
line direction usually is not uniform (Fig.1.18], Some low fold areas (holes] 
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are created, which must be filled in during acquisition by shooting more 
lines at appropriate locations. 

Seismic Integrated Positing System fSIPSl 

The SIPS is an acoustic positioning system for submerged parts of the 
seismic configuration. For multi-streamer and multi-source, it is very 
essential to know the dynamic separations between cables, between sources 
and between source & cable for proper binning of the data. The SIPS fulfils 
this need. Basically it measures the ranges between the deployed acoustic 
units using the principle that sound travels through water at a predictable 
speed. By measuring the time between transmission and arrival of the 
signal, the distance between two units can be calculated. The sound 
generated by SIPS units is short acoustic pulses of a constant high frequency. 
This frequency is different from the frequency of the airgun-generated signal 
or the signal received by the hydrophones in order to prevent interference 
in the seismic data. 

The SIPS makes use of transceivers for the range determination that consist 
of a transmitter, a receiver and equipment to measure the travel time. There 
are two sub-systems, which make up the total SIPS system. They are the 
following: 

The Cross Streamer Ranging System (XSRS] 

The Head And Gun Position System (HGPS] 

The XSRS units are attached to the streamer and are able to communicate 
with the ship borne equipment via a digital data link through the coils in 
streamer and XSRS units. Since the acoustic transducer is placed in the 
vertical axis, an omni-directional pattern is achieved which means that 
ranging can take place up and down the streamer as well as between 
streamers. 

The HGPS transceivers are mounted on the gun arrays and the hull. They are 
more rugged to sustain the jerk and vibrations of ship and airgun 
respectively. They are powered by ship-borne sources and they 
communicate with the onboard equipment through a pair of wires. 

The SIPS system is further divided in the following three separate acoustic 
nets. 
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Front net, Middle net and Tail net 

Net lay outs for six streamers are given below in Fig. 1.19 
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Fig. 1.19: Net layout diagram for six streamers 

The front net determines ranges between the transceivers mounted on the 
hull, airgun arrays and the two front XSRS units of each streamer. The HGPS 
units are positioned relative to the vessel through their range measurements 
by which the position of the XSRS units can be determined. 

The middle net consists of two XSRS units located at the middle of each 
streamer. Since there is no absolute positioning system near this location, 
only the mid separation of the streamer can be determined from this data. 
Furthermore the data can be compared with the various compass readings 
near the XSRS units. 

Just like the front net, the tail net also contains absolute reference points, 
which are the tailbuoys. The tail XSRS unit provides the signal that triggers 
the transponder at the tail-buoy. 

Ranges to and from the head network gun transceivers sould be measured 
either before the guns fire or as long after as practicable to avoid the 
mechanical shock wave generated. 

The following are the typical problems which sometimes occur: 

False echoes : False echoes arise when there is a good reflector in shallow 
water. The transceivers will then receive two signals soon after each other, 
the direct and the reflected signal 
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Velocity-distortion: Velocity distortion occurs when measurements are 
executed immediately behind the vessel and gun strings. At these locations 
the water contains a relative large amount of air bubbles that reduces the 
speed of sound. The area affected by the propellers and guns are small since 
the bubbles soon emerge at the water surface. 

For every shot-point one cycle of range measurement is executed, which 
consists of several active events. The three nets operate independently of 
each other and data can be collected simultaneously in one event. The 
controller of the SIPS is informed when a shot-point is reached and after a 
brief delay, it commands the XSRS and HGPS controllers to start their range 
measurements. After completing the cycle, the controllers request all the 
ranges measured by the transceivers. The relative positions of the 
transceivers are thus known with a precision of 0.05 m. by tying in these 
positions with the absolute positions of the vessel and tail buoys, the true 
positions can be known. 
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CHAPTER 2 


GEOMETRICAL PARAMETERS 


2.0 INTRODUCTION 

G eometrical parameters mainly consist of far offset, near offset, array 
length and group interval, which affect the quality of the data in 
different ways. Far-offset should be large enough for better velocity analysis 
and multiple attenuations but it should be within a reasonable limit, so that 
NMO stretch should not be more. Near offset should be such that signal does 
not get lost amidst shot generated noise. Array length should be as small as 
possible depending upon the noise characteristics, so that shorter 
wavelength component of signal, required for better resolution, may not fall 
into the reject band of the array. Group interval is nothing but sampling 
interval in space domain. Therefore sampling theorem must be satisfied for 
the shortest apparent wavelength component of the signal, which we want 
to record as per geological objectives. There are other factors like, 
multiplicity, line length, line orientation & line spacing, which will also be 
discussed in this Chapter. 


2.1 FAR OFFSET DISTANCE 

In reflection seismology, we record the two way time (t] of reflected events 
as a function of offset distance (x). As the offset distance increases, two way 
time also increases. For a simple case of plane reflector in which overlying 
medium has uniform velocity V, the relation between two way time [t] and 
offset distance (X] will be given by (Fig. 2.1]. 

V 2 t 2 = V 2 t 0 2 + X 2 

t 2 = X 2 /V 2 + t 0 2 (2.1) 

t 2 / to 2 - x 2 /vV = 


1 
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Where, t 0 = two way time for the normal reflection. 

This is the equation of hyperbola. 

The time difference (t - t 0 ], the excess time between slant and normal two 
ways reflection time is called normal move out (NMO) 



Fig.2.1: Ray path Geometry 
In the above equation (2.1] if we put, 


AT =t-t 0 

We get, from equation [2.1] 

(AT + t 0 ) 2 = t 0 2 + X 2 / V 2 
Or (AT) 2 + 2 AT.t 0 = X 2 / V 2 

For small offset distance, AT being small, if we neglect (AT) 2 in the above 


equation, we find 


AT= X 2 /2t 0 V 2 

(2.2) 

Or, X 2 = 2t 0 V 2 AT 

(2.3) 

Or, V 2 = X 2 / 2t 0 AT 

(2.4) 
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While processing the data, normal move-out correction is applied to the 
signal, i.e. AT time is subtracted from the two way time, to convert the two 
way slant time into two way normal reflection time. 

The signal is not a spike, but it is a pulse of finite duration. Since Normal 
Move out 'AT' is inversely proportional to time 'to', both end of the pulse will 
not have the same correction, early part of the signal (corresponding to 
lower to] will have more correction than the later part (corresponding to 
larger to]. Therefore width of the signal will increase after NMO correction as 
shown in Fig. 2.3. Thus higher frequency signal changes into slightly lower 
frequency after NMO correction. This effect is called NMO stretch. As we will 
see, the larger the offset the higher will be the stretch in signal and thereby 
we will get lesser resolution. Thus to have better resolution of beds, we must 
not exceed the far offset beyond certain maximum value. 

According to equation (2.2], for the determination of AT, we must know the 
velocity V. Therefore before NMO correction, velocity analysis is carried out. 
From equation (2.1] X 2 vs. t 2 curves is a straight line, intercept of which on 
the time axis gives 'to 2 ' and inverse of the slope gives the value of V 2 . But in 
real situations, there is not only one interface as shown in Fig.2.1 but also 
large number of interfaces separating beds of different acoustic impedances 
and at every interface reflection and refraction both take place and hence 
ray path geometry is not a straight line for all the interfaces. Consequently 
time distance relationship (hodograph] will not be true hyperbolic, but it 
will only approximate to a hyperbola. Therefore t 2 vs. X 2 plot will not be a 
straight line, but a curve slightly concave towards the X 2 -axis. The square 
root of inverse of the slope of the best-fit straight line passing through this 
curve is called stacking velocity 'Vs'. This velocity is used for NMO correction 
prior to stacking. In practice stacking velocity search is carried out within a 
time gate generally equal to the dominant period of the reflection. Several 
velocities are tried in the gate for NMO correction one by one and NMO 
corrected samples are summed within the gate. The velocity producing 
maximum sum in a gate is the stacking velocity for that gate. If NMO of the 
events are such that the hyperbolic events are well within the gate, a small 
variation in the stacking velocity will not alter the value of summation and 
hence there will be loss in velocity resolution. Therefore to have better 
velocity resolution hyperbolic event must fall outside the time gate used for 
velocity analysis. After NMO correction all the events along one hyperbola, 
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will align along a straight line perpendicular to time axis and when stacked 
(or summed] all traces will be added in phase, producing maximum 
amplitude. On the other hand, multiple events travel in upper medium with 
lower velocity (in normal cases] NMO of the multiple event will be larger 
than primary event at the same offset and therefore stacking velocity for 
NMO correction required to align these primary events perpendicular to 
time axis, will not be able to align the multiple events and when summed 
multiple event will not be added in phase and hence will be of lesser 
amplitude. This is the principle behind multiple attenuations. If difference 
between NMO of P & M (P-Primary, M-Multiple] is very small, we will not be 
able to differentiate primary and multiple events, separately, instead we will 
get an interference pattern of these two events producing distortion in wave 
shape and attenuation of multiple will not be possible. For identification & 
attenuation of multiples, the difference of NMO between P & M must be at 
least equal to the period of the reflected event. 

Thus the requirement of good velocity analysis and multiple attenuation 
emphasizes the need of larger offset distance, but to avoid larger NMO 
stretch, shorter far offset distance is desirable. Therefore, far-offset distance 
must be decided carefully to get optimum result. 

2.1.1 FAR OFFSET DISTANCE FOR VELOCITY RESOLUTION 

It was discussed above, that velocity search is carried out inside a time gate. 
For better velocity resolution, hyperbolic event must fall at least just outside 
the gate at farthest trace. In other words, NMO of the trace at farthest offset 
must be at least equal to or greater than the time gate of velocity search. 
Time gate of velocity search is usually taken equal to the dominant period of 
reflection event. Thus if we put At= 1.5 T, in equation 2.2, where T is the 
dominant period of reflection, we shall get 


X min 2 = 3 V 2 t 0 T 
X mi n = (6ZVT) 1/2 
Xmin = V (3 to T) 1/2 


(2.5i) 

(2.5ii) 


Where Z = Vto /2, the depth of target horizon. Equation 2.5 (i] & (ii] give the 
minimum value of far - offset distance for better velocity resolution. 



41 


Example 2.1 

Calculate the minimum far offset distance to map a reflector where t 0 is 1.6 
sec., velocity is 2500 m/sec., and dominant frequency of reflection is 20 Hz 
(i.e. T = 50ms] 

Solution: 

From equation (2.5ii) 

X = 2500 (3 x1.6 x0.05)1/2 = 1225 m 

2.1.2 FAR OFFSET DISTANCE FOR MULTIPLE ATTENUATION 

For effective attenuation of multiple, the differential move out (difference of 
move outs) of primary and multiple reflections should be at least equal to or 
greater than the period of the multiple event at far offset distance, 
preferably 1.5 times the period of multiple or time gate of coherency 
measurement to be used for velocity analysis. Normal Move-Outs for 
primary and first order multiple reaching around the same time 'to is given 

by 

(AT) P = (t 0 2 + X 2 / V 2 2 ) 1/2 - t 0 

(AT) M = 2 [{(t 0 /2) 2 + (X/2) 2 / Vi 2 } 1/2 -t 0 /2] 

(At) m = (t 0 2 + X 2 / Vi 2 ) 1/2 -t 0 



Fig.2.2: Variation of NMO with offset 
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Where, 'to' is the two way normal reflection time for primary event and Vi & 
V 2 are the velocities corresponding to multiple & primary reflections. A 
typical curve for NMO of primary & multiple reflections is shown in Fig. 2.2 

Differential Move-out is given by 

(AT] d = (AT] M - (AT] P 

= (t 0 2 + X 2 / V x 2 ) 1/2 - (t 0 2 + X 2 / V 2 2 ) 172 

= t 0 (1 + X 2 / 2t 0 V! 2 ) - t 0 (1 + X 2 / 2t 0 V 2 2 ) 

= [X 2 /2t 0 ] [1/ V x 2 - 1/ V 2 2 ] 

Values of (AT]p, (AT] m and (AT]d have been calculated and shown in Table- 
2.1 


Table 2.1 


X 

(AT] P 

(AT] m 

(AT]d 

500 

0.011 

0.013 

0.002 

1000 

0.043 

0.052 

0.009 

1500 

0.097 

0.117 

0.020 

2000 

0.172 

0.208 

0.036 

2500 

0.269 

0.326 

0.057 

3000 

0.387 

0.469 

0.082 


If we put (AT]d = 1.5 T, we shall get 

X 2 = 3 t 0 T/ [1/ V x 2 - 1/ V 2 2 ] (2.5 iii) 

Where T = Time period of the signal of smallest frequency of interest. 

Example 2.2 

To calculate far offset distance, if Vi = 2000 m/sec., V 2 = 2200 m/sec. 
to = 2.4 sec. and T = 50 ms. 

Solution: 

From relation, equation 2.5 (iii], we get X = 2880 m. 

The same result may also be obtained from the primary & multiple move out 
curves in Fig. 2.2. 
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2.1.3 FAR OFFSET DISTANCE CONSIDERING NMO STRETCH 

As discussed earlier, stretching in signal is produced as a result of unequal 
NMO correction in early and later parts of the wavelet (Fig.2.3). Amount of 
stretching depends on the rate of change of NMO with respect to time. 

NMO is defined as 

AT = t - t 0 

Where, t = two way slant reflection time & t o= two way normal reflection 
time. 

Therefore, 

d/dt (AT) = 1 - (dto/dt) (2.5 iv) 

d/dt(AT) is always negative, because as't' increases AT decreases . 

If d/dt (AT) would have been zero or signal would have been ideal spike 
there would not have been the problem of signal stretch. Therefore, NMO 
stretch factor is defined as : 

S = - d/dt(AT) 

= (dto/dt - 1) 

For a single non-dipping reflector, since 
t = (t 0 2 + X 2 / V 2 ) 1/2 
dt = t 0 dto/ (to 2 + X 2 / V 2 ) 172 
dto/dt = (t 0 2 + X 2 / V 2 ) 1/2 / t 0 

dt 0 / dt = Sec 9 

Where, 9 is the angle between reflected wave and normal to the reflector, 
(i.e. angle of reflection) 

Hence. S = (Sec 9 -1) 

In terms of percentage, S = (Sec 9 -1) x 199% 

Stretch may also be expressed as: 

S = [(t 0 2 + X 2 / V 2 ) 172 / t 0 ] - 1 
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Fig. 2.3: NMO Stretch 
S = ( 1+ X 2 /V 2 t 0 2 ) 172 -1 

S = ( 1 + X 2 /4Z 2 ) 1/2 - 1 (2.5 v) 

Where, Vt 0 = 2Z 

Variation of 'S' with respect to (X/Z] has been shown in Fig 2.4. We find 
from above, that 

X 2 /4Z 2 = (S+l) 2 - 1 

X = 2Z ( S 2 + 2S) 1/2 

If we put, S = 0.10 (i.e. 10% stretch], in the above equation, we get 

X = 0.92 (Z) (2.5 vi) 

Thus as a rule of thumb far offset distance should not be kept more than the 
depth of target horizon to have less than 10% stretch. 

Table 2.2 


X/Z 

S (%] 

0.5 

3 

1 

12 

1.5 

25 

2 

41 

2.5 

60 

3 

80 
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Table-2.2 and Fig.2.4 show the variation of stretch with [X/Z] 

Example 2.3 

Considering depth of shallowest reflector Z = 500 m, calculate stretch [a] 
when X = 1000 m. (b] when X = 800 m. 

Solution: 

(b] when X = 1000 m. 

stretch = 41% 

(b] when X= 800 m. 
stretch = 28% 


2.2 NEAR-OFFSET DISTANCE 

Near Offset Distance is decided on the basis of the noise interference with 
near traces. It is better to select the minimum near offset to preserve the 
proper foldage for shallow reflections. Near trace offset should not be so 
large that shallow reflections are muted due to NMO stretching at processing 
centre. If signal to shot generated noise ratio is very less and signal is not 
recoverable from the noise due to distortion produced by the recording 
system, we may choose larger near trace-offset otherwise, near trace offset 
should be as short as possible. 

Other reasons for keeping minimum offset close to zero are to get better 
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velocity and timing control. It is good to have a near-zero offset data for 
processing of first arrival information near the source for static & dynamic 
correction. When offset is large, these events are not recorded. 

Considering the example 2.3, if stretching of more than 30% is muted in 
processing centre, only the channels situated up to 800 m from SP will be 
retained for the signal from the reflector at depth of 500 m. If group interval 
is 50 m and near-trace offset is also 50 m, then only 16 channels are 
recording the useful signal and only 8 fold data is being recorded. If near 
offset is kept 300 m, then only 10 channels between X = 300m and X = 800m 
is effectively recording the signal and only 5 fold data can be obtained for a 
reflector at Z = 500 m. Thus as near offset is increased, we get data of lesser 
foldage at shallower level. 


2.3 ARRAY LENGTH AND EFFECT ON RESOLUTION 

Desired shortest wavelength component of signal should always lie within 
the pass band (attenuation not more than 3db] of the array for better 
resolution and attenuation of noise must be more than 20 db (reject band] 
for better signal to noise ratio. For any linear array of equal weight and 
equal spacing, position of pass & reject bands depend on the effective array 
length (L = N.d], where N is number of elements and d is element spacing. 
Let the desired shortest apparent wavelength of signal at far offset be Xsx & 
longest wavelength of noise be Xn. Then for better resolution and good signal 
to noise ratio, Xsx must not be attenuated more than 3 dB & Xn must be 
attenuated at least by 20 dB by the array. For any linear array of equal 
weight & equal spacing (See Chap. 4], 3 dB attenuation approximately occurs 
at L/Xsx = 0.44 & 20 dB attenuation approximately occurs at L/Xn= 0.91. 

Thus we must have (See Fig. 2.5] 

L < 0.44 X sx 

L > 0.91 X N 


( 2 . 6 ) 

(2.7) 
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Fig.2.5: Respose curve for 6 elements uniform arrray 



Since, apparent velocity Vx =V / SinO, (Fig. 2.6], therefore 
Ax = A/ Sin 9 and 
SinO = 7,s/7.sx 

Putting the limiting value of Ax = L/0.44 form equation (2.6) 

SinO = 0.44 A /L (2.8) 

Putting the limiting value of L = 0.91 In from equation (2.7) 

SinO = 0.48A An 


(2.9) 
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Since for plane reflector, tan 0 = X/2Z, therefore, Sin 0 = 1/ (1+4Z 2 /X 2 ) V 2 
and hence from equation (2.8] by putting the value of Sin0 

(1+4Z 2 /X 2 ) 1/2 = L/ 0.44X S , 

L = 0.44X S (1+4Z 2 /X 2 ) 1/2 (2.10) 

Which gives 

X/Z = 2/[(2.27 L/Xs) 2 -l] 1/2 (2.11) 

And from equation [2.9], putting Sin 0 = 1/ (1+4Z 2 /X 2 ) : /2 
(1+4Z 2 /X 2 ) 1/2 = X N / 0.48X S , 

This gives 

X/Z = 2/ [{2.08X N /X S } 2 - 1 ] 1/2 ( 2 . 12 ) 

Value of L for different value of X/Z and value of X/Z for different values of 
the ratio Xn/Xs have been calculated using equations [2.10] & [2.12] and 
shown in Table- 2.3 & 2.4 and corresponding plots are shown in Fig. 2.7 & 
2.8 respectively. 


Table 2.3 


X/Z 

Effective Array Length L [m] 

1 

2 

3 

4 

5 

Xs=20 m 

Xs=30 m 

Xs=40 m 

Xs=50 m 

Xs=80 m 

0.5 

37 

54 

74 

91 

144 

0.6 

33 

47 

66 

80 

127 

0.7 

28 

41 

57 

70 

110 

0.8 

24 

35 

48 

59 

94 

1.0 

20 

29 

40 

49 

78 

1.3 

16 

23 

32 

40 

63 

2.0 

13 

18 

25 

31 

49 

4.0 

10 

15 

20 

26 

39 
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Table 2.4 


An/Xs 

X/Z 

0.6 

2.7 

0.7 

1.9 

0.8 

1.5 

0.9 

1.3 

1.0 

1.1 

1.1 

1.0 

1.2 

0.9 

1.3 

0.8 

1.4 

0.7 
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Fig.2.8: relation between A/A and X/Z 


2.4 EFFECT OF GROUP INTERVAL ON RESOLUTION 

Group Interval is the spatial sampling interval in receiver domain. Therefore 
condition imposed by sampling theorem must be satisfied for true recovery 
of the signal from sampled values. Sampling theorem states that, there must 
be minimum two samples per cycle of the wave to be recovered. If 0 is the 
emergent angle, then apparent wavelength corresponding to shortest 
wavelength of the signal that is required for recording is 

Ax = A/ Sin 9 [Fig. 2.6] 

Where, A = Shortest wavelength of signal desired to fulfil the objective of 
the survey. 

For plane reflector, tan 9 = X/2Z, Where, X = far offset distance, Z = depth of 
the target horizon, then 

Sin 9 = 1/ (1+4Z 2 /X 2 ) 1/2 and Ax = A (1+4Z 2 /X 2 ) 1/2 

Therefore, maximum Group Interval that can be possible as per sampling 
theorem will be 

AX = Ax / 2 and therefore, 

AX = (0.5 A) (1+4Z 2 /X 2 ) 1/2 


(2.13) 
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Variation of Group Interval (AX) with respect to (X/Z) for different values of X s 
has been calculated in the following Table-2.5 and shown in Fig. 2.9. 

Table 2.5 


X/Z 

Max. Group interval (m) 

1 

2 

3 

4 

5 

Xs=20 m 

Xs=30 m 

Xs=40 m 

Xs=50 m 

Xs=80 m 

2.0 

14 

21 

28 

35 

56 

1.67 

16 

24 

32 

40 

64 

1.25 

19 

28 

38 

47 

76 

1.0 

22 

33 

44 

55 

88 

0.83 

26 

39 

52 

65 

104 

0.71 

30 

45 

60 

75 

120 

0.63 

34 

51 

68 

85 

136 

0.56 

37 

55 

74 

92 

148 

0.5 

41 

61 

82 

102 

164 



Fig. 2.9: Plot of maximum group interval vs. X/Z 
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2.5 EFFECT OF DIP ON SPREAD GEOMETRY 

In the earlier consideration, it has been assumed that the reflector is plane 
and there is no dip. Now we shall consider effect of dipping reflector on [i] 
emergent angle & incident angle [ii] move-outs (hi] far offset distance (iv) 
array length and [v] group interval. 

2.5.1 EFFECT OF DIP ON EMERGENT AND INCIDENT ANGLE 

(A) Dipping Reflector: Down-Dip Energy Transmission 
In Fig.2.10 

0 = Emergent angle for the ray reaching the detector at far offset [X] 
a = Incident angle for the ray reaching the detector at far offset [X] 

® = Angle of dip 
0S’= 2Z Cos® 

OS = 2Z Sin ® 

Tan 0 = (X + 2Z Sin ®]/ (2Z Cos®] 

AX/At = V/ Sin (a +® ] 

Tan 0 = (X+2Z Sin ®]/(2Z Cos®] (2.14] 

0 = a+ ® (2.15) 



Fig, 2.10: Ray Geometry in Up-Dip shooting 
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(B) Dipping Reflector: Up-Dip Energy Transmission 
In Fig.2.11 

0 = Emergent angle for the ray at far offset (X] 
a = Incident angle for the ray at far offset (X] 

® = Angle of dip 
0S’= 2Z Cos® 

OS = 2Z Sin ® 

Tan 0 = (X - 2Z Sin ®]/ (2Z Cos®] 

AX/At = V/ Sin (a -®] 

tan 0 = ( X-2Z Sin ®]/(2Z Cos ® ] (2.16) 



Fig. 2.11: Ray Geometry in Up-Dip shooting 
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Variation of emergent angle ' 0 ' with (X/Z] for different value of dip angle 
' ® ' has been shown in Fig. 2.12 and also in the following Table 2.6 

Table 2.6 



0 

0 

0° 

0 

0 

X/Z 

Up-dip 

No dip 

Down-dip 


O 

O 

t-H 

II 

e 

o 

LD 

II 

e 

O 

O 

II 

e 

o 

LD 

II 

e 

® =10° 

0.5 

4 

9 

14 

19 

23 

0.6 

7 

12 

17 

21 

26 

0.8 

13 

17 

22 

26 

30 

1 

18 

23 

27 

31 

34 

1.2 

23 

27 

31 

35 

38 

1.4 

28 

32 

35 

38 

42 

1.5 

30 

34 

37 

40 

43 

1.6 

32 

36 

39 

42 

45 

1.8 

36 

39 

42 

45 

47 

2 

40 

43 

45 

48 

50 



Fig.2.12: Variation of emergent angle with X/Z 
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EXAMPLE2.4 


Let us consider one practical example for a reflector at depth. Z = 2500 m. & 
offset X = 2500 m. for calculation of emergent and incident angle for (a] non¬ 
dipping reflector (® =0] and (b] dipping reflector (® = 10°) for up-dip case 
& (c] dipping reflector (® = 10° ] for down-dip case. 

Solution: 


a = 26 ° 


(a) Non-dipping reflector (®=0°) 

Emergent angle 0 = 26° 
Incident angle 

(b) Down dip (® = 10°) 

Emergent angle 
Incident angle 

(c) Up dip (®= 10°) 

Emergent angle 0 =18° 
Incident angle a =28° 


[from eqn.2.14] 
[from eqn.2.15] 


0 = 34° [from eqn.2.16] 
a=24° [from eqn.2.17] 


[from eqn.2.16] 
[from eqn.2.17] 


Thus we see that emergent angle is less in the case of up-dip and more in the 
case of down dip shooting, which implies that apparent wave length of signal 
will be more in the case of up-dip and hence we may dare to have longer far 
offset and larger group interval in this case than in the case of down dip 
shooting without having undesirable effect of array on the signal. 


2.5.1 EFFECT OF DIP ON MOVE OUT. 


(A) Dipping Reflector: Down Dip Case 

Move out time = t -1 o, where 
t 2 = (2Z Cos ® / V) 2 +{ (X + 2Z Sin ®) / V} 2 
t 0 =2[{Z + (XSin®)/2}/V] 
t & t o have usual meaning 
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(B) DIPPING REFLECTOR: Up-Dip CASE 

Move out time = t -1 o, where 
t 2 = (2Z Cos O / V) 2 +{ (X - 2Z Sin O) / V} 2 
t 0 = 2 [{ Z - (X Sin O) / 2 }/V ] 
t & t o have usual meaning 

Example2.5 

Now let us consider one practical example for a reflector at depth. Z = 2500 
m. & offset X = 2500 m. for (a] non- dipping reflector (® =0] and (b] dipping 
reflector (® = 10° ] for up-dip & (c] dipping reflector (® = 10° ] down-dip 
case to calculate move out time. 

Solution: 

(a] Non- dipping reflector (® =0] 

t = 2.236 Sec., t o = 2.000 Sec., 

AT = 0.236 Sec. 

(b] Dipping reflector (® = 10° ] for down-dip case 

for above example- 
t = 2.386 Sec. 
t o = 2.174 Sec. 

AT = 0.212 Sec. 

(c] Dipping reflector (® = 10° ] for up-dip case 

for above example- 
t = 2.075 Sec. 
t o = 1.826 Sec. 


AT = 0.249 Sec 
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2.5.2 EFFECT OF DIP ON FAR OFFSET DISTANCE 

(A) Dipping Reflector: Down-Dip Energy Transmission 
In Fig.2.10 

X=SRi=OR r OS 

And since, ORi=2Z Cos® tanG and OS= 2Z Sin® 

X = 2Zcos®tanO - 2Z Sin® (2.18) 

(B) Dipping Reflector:UP-Dip Energy Transmission 
In Fig.2.11 

X=SRi=ORi+ OS 

X = 2Zcos®tanO + 2Z Sin® (2.19) 

Thus for the same value of emergent angle 0 (i.e. for same apparent 
wavelength and in turn for the same array length] we can have more far 
offset in the case of Up-dip than in down dip case. 

Example 2.6 

To calculate far offset distance, if 0= 30° as calculated using equation Sin 0 = 
0.44 X s / L, be the maximum limit that can be had on the basis of array 
length and shortest wavelength of signal for a reflector where, Z= 2500 m 
and 0=10°. 

Solution: 

Considering the previous example: (Z = 2500 m, ® = 10 °) 

(a) Non- dipping reflector (® =0) 

X = 2886 m 

(b) Dipping reflector (® = 10° ) for down-dip case (equation 2.18) 

X= 1975 m 

(c) Dipping reflector (® = 10° ) for up-dip case (equation.2.19) 
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X = 3710 m 

2.5.3 EFFECT OF DIP ON ARRAY LENGTH 

As discussed in section (2.3], effective Array length is given by 
L = 0.44 A, s / Sin 0 

(A) Non-dip case 

Substituting the value of Sin 0 from equation 2.10, effective array length is 
L n = 0.44 A, s ( 1 + 4Z 2 / X 2 ) 1/2 (2.20) 

Where, Ln = Effective array length in non- dipping case 

(B) Down-dip case 

Substituting the value of Sin 0 from equation 2.14 for down-dip case, 
effective array length is 

L d = 0.44A. s [ 1 + {4Z 2 Cos 2 O } / { X + 2Z Sin ® } 2 ] 1/2 (2.21) 

Where, L D = Effective array length in down- dipping case 

(C) Up-dip case 

Substituting the value of Sin 0 from equation 2.16, for up-dip case effective 
array length is 

Lu = 0.44iA s [ 1 + {4Z 2 Cos 2 O } / { X - 2Z Sin O } 2 ] 1/2 (2.22) 

Where, Ly = Effective array length in up- dipping case 

Thus, we see that, L u > Ln > Ld 

Values of (L /A.s) for different angle of dip <j> has been calculated in Table- 
2.7 and plot is shown in Fig. 2.13. 
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Table 2.7 


X/Z 

L/As 

(j> = 0 o 
(NO-DIP) 

(j>= 10 o 
(UP-DIP) 

(() =10° 
(DOWN-DIP) 

0.5 

1.8 

5.7 

1.1 

0.6 

1.5 

3.5 

1.0 

0.8 

1.2 

2.0 

0.9 

1.0 

1.0 

1.4 

0.8 

1.2 

0.9 

1.1 

0.7 

1.5 

0.7 

0.9 

0.6 



Fig.2.13: Variation ofL/As with X/Z 
2.5.4 EFFECT OF DIP ON GROUP INTERVAL 

As discussed in section (2.4], group interval is [AX] =0.5 A s / Sin 0 

(A) Non-dip case 

From eqn.2.13, for non-dip case, group interval is, 

(AX) N = 0.5 A s (1+4Z 2 /X 2 ) 1/2 


(2.23) 
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(B) Down-dip case 

Substituting the value of Sin 0 from eqn.2.14, for down-dip case, 

(AX) D = 0.5 X s [1+ {4Z 2 Cos 2 ®}/ {X+2Z Sin®} 2 ] 1/2 (2.24) 

(hi) Up-dip Case 

Substituting the value of SinO from eqn.2.16, for up-dip case, group interval 
is 

(AX) u = 0.5 X s [1+ {4Z 2 Cos 2 ®} / {X-2Z Sin®} 2 ] 1/2 (2.25) 


Thus we see that 
(AX) u > (AX) N >(AX) D 

Variation of (AX /X s ) with angle of dip (]) has been calculated in Table- 2.8 
and plot is shown in Fig. 2.14. 

Table 2.8 


x/z 

AX As 

<J>= 10° 
(DOWN-DIP) 

<]>= 0 o 

(NO-DIP) 

()>= 10 o 

(UP-DIP) 

0.5 

1.31 

2.0 

6.4 

0.6 

1.1 

1.7 

4.0 

0.8 

1.0 

1.4 

2.3 

1.0 

0.9 

1.1 

1.6 

1.2 

0.8 

1.0 

1.3 

1.5 

0.7 

0.8 

1.0 


Thus it is seen that in the case of up-dip shooting we may go for slightly 
higher values of group interval than in the case of down dip shooting for the 
same far offset & resolution. Therefore, generally up dip shooting is 
preferred. 
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2.6 FOLD AGE REQUIREMENT 

Multiplicity or foldage is the number of times each subsurface reflection 
point is sampled through different ray-path geometry in CDP profiles. 
Amount of multiplicity required can be determined only by comparing the 
quality of sections processed from the same data with different multiplicity. 
In general 600% multiplicity is the minimum requirement for proper 
functioning of processing routines, which require multiple traces at each 
CDP. Doubling the multiplicity gives a theoretical improvement of 6 db in 
signal to random noise ratio and also better rejection of the other noises. In 
most areas multiplicity greater than 2400% give only marginal 
improvements. When targets of interest are both shallow and deep, muting 
reduces the multiplicity at shallow target; hence high multiplicity in this case 
is most useful so that minimum needed multiplicity at shallow target is 
retained even after muting. Once the number of recording channels is set by 
the consideration of near-offset, far-offset and group interval the amount of 
multiplicity sets the shot point spacing in terms of geophone group interval. 

Multiplicity = (1/2) x (No. Of channels x Group interval) / (Shot point 
Interval) 

Provided, shot- interval >_Group - interval 
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2.7 LINE-LENGTH AND RECORD-LENGTH 

Event at point A migrated to point Am (Fig.2.15] 
(j) = Apparent Dip 
4>m= True Dip after Migration 
Sin 4) M = Tan 4> 

X Xm = Migration Aperture 



For effective migration, a line should extend beyond the area of interest by a 
distance at least equal to the migration aperture (XX m = Z tan cj)] (Fig. 2.15] 
for the deepest target at both ends. Migration aperture is the horizontal 
distance between events at its apparent & true positions. 

Example 2.7 

To calculate Migration Aperture, if depth of the target horizon is Z = 3000km 
and angle of dip <J) = 40 0 (apparent dip]. 

Solution: 

Migration aperture, 

XX m = Z tan (]) = 3000 X (Tan 40 °] m = 2520 m (approx.] 
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Therefore line length should be increased by 2520 m at both the ends. 

From Fig.2.15 we see that the event A at time To has shifted to A' at time Tom. 
Therefore if zone of interest is at time Tom record length must be kept To for 
proper migration of the event where 

To = Tom Sec 4 >m. 

Also the line should have full foldage at both the ends. To fulfill this 
requirement we would have to take extra shots either at both ends or at one 
end depending on the type of shooting, split-spread or end-on. Whatever the 
shooting technique, last CDP of first record corresponds to full fold CDP at 
the beginning of the line and first CDP of the last record corresponds to full 
fold CDP at the end of the line. BC is the subsurface length of the line having 
full foldage and AD is the actual length of the line on the ground as shown in 
Figs. 2.16, 2.17 and 2.18. 
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2.8 ORIENTATION OF SEISMIC LINES 

The more reliable interpretation can be made of a dip line than of a strike 
line or one oblique to the dip and also conventional migration depends on 
the assumption that the line is in dip direction. Therefore, line orientation 
should be decided on the basis of the geological structure to be mapped. The 
mapping along strike lines, perform two functions. Firstly, they confirm the 
correlation between dip lines and secondly, they show that the dip lines are 
indeed dip lines. In the area where the direction of dip varies with depth or 
where strikes changes rapidly, the distinction between dip and strike lines 
may not be clear. All directions are equally important in this case and this 
leads to a requirement of three-dimensional recording. Dip lines should be 
close enough to detect the smallest prospective feature. Strike lines can be 
spaced at a larger distance than dip lines. 


2.9 SPACING OF SEISMIC LINES 

For a reconnaissance survey the distance between the two lines may be kept 
such a large that a general sub-surface feature of interest cannot be missed. 
For detailed surveys, the dip lines should be kept close enough so that 
variation in strike direction can be sampled adequately, where the strike 
changes rapidly. At the time of contouring if there is trouble in deciding that 
contour may go between the lines, it indicates that the lines are too far apart. 


2.10 PRACTICAL EXAMPLE FOR SPREAD DESIGN 
EXAMPLE2.8 

To design the spread geometry for an area, where 
Minimum Depth of shallowest target = 2500 m (main reflector) 

Depth of deepest target = 7000 m 
Maximum Depth of main reflector = 3000 m 

Average velocity of the shallowest target at minimum depth =2400 m/sec 

Average velocity of the deepest target =3000 m/sec 

Average velocity of the main reflector at maximum Depth =2500 m/sec 
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Main objective of the survey is to map sand bars of thickness 10 m to 20 m at 
main reflector level. 

Solution: 

Two way normal reflection time at shallowest level 
= 2 (2500/2400] =2.080 sec. 

Two way normal reflection time at deepest level = 2 (7000 / 3000] 

= 4.670 sec. 

Two way normal reflection times at Maximum Depth of the main reflector 
= 2 (3000 / 2500] = 2.4sec. 

Max. Dip at main reflector = 3 0 

(1) Shortest Wavelength of Interest (A, s) at Main Reflector Level 

As per objective of the survey, minimum thickness of the object to be 
resolved = 10 m. Considering that resolvable limit is A s /4, therefore, 

A s /4 = 10m. A s = 40 m. 

This is the required shortest wavelength of signal to fulfil the geological 
objective. If V=2400 m/sec, highest frequency of interest f h = 2400 / 40 = 60 
Hz 

(2) Far offset distance 

(a] Requirement for Better Velocity Analysis At Deepest 
Reflector Level 

As per objective, our deepest target is at 7000 m; therefore, at least 
velocity analysis must be reliable up to this level. Hence putting Z = 
7000 m, V=3000m/sec, T = 50 ms. (Time period of event at this level 
assuming 20 Hz signal at this level] 

X min = ( 6ZVT)% = 2500 m. [Eqn.2.5(i]] 

(b] Requirement for Multiple Attenuation at Deepest Part of 
Main Reflector Level 
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Suppose that there is multiple activity in the area at main reflector 
level. Therefore putting Z = 3000 m (deepest part], V 2 = 2500 m/sec. 
(assumed], to = 2.4 sec. and Vi = 2300 m /sec. (assumed] 

X 2 = [ 3 t 0 T/ (1/ Vi 2 - 1/ V 2 2 } ] % = 2950 m [Eqn.2.5(iii]] 

(c] Stretch Effect at Shallowest Level of Main Reflector 

Stretch S = (1 + X 2 /4Z 2 ] *4 - 1 [Eqn.2.5(V]] 

X = 2950m (as calculated at 'b' above], Z = 2500 m, 

X/Z = 1.18 

S = 0.16 ( = 16%] 

Stretch is more than 10%, and therefore some of the farther traces 
beyond X = 2500 m ( 15 traces, if group interval 30 m is assumed ] 
may be muted. But these traces may provide valuable guide for 
velocity analysis & multiple attenuation. 

(3) Effective Array Length ( L = Nd) 

Now, X s = 40 m (as already calculated] 

And X/Z = 1.18 (as already calculated] 

L= 0.44X s (1 + 4Z 2 /X 2 ) % = 35 m [From eq.2.10] 

Thus we have to choose only those geophone pattern which have 
Nd = 35m (this is for uniform array only, tapered array may have 
different effective length]. Value of N & d may be changed according 
to characteristics of surface noise. 

(4) Maximum Group Interval ( AX ) 

Since X s= 40 m and X/Z = 1.18 

AX = 0.5 X s (1 + 4Z 2 /X 2 ] % = 39 [From eq. 2.13] 

Since angle of dip is small, we have not considered the dip effect. 
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CHAPTER-3 


SOURCE PARAMETERS 


3.0 INTRODUCTION 

P roper source signature is most essential factor to fulfill the geological 
objective of the survey. Because of the two important requirement of 
good data quality (i) the resolution and (ii) the signal to noise ratio are 
mainly dependent on the amplitude and frequency of the source wavelet. 
Spike is the ideal source signature but it is impossible to achieve in practice. 
Any type of source whether Explosive or Vibrator or Air gun generates 
different kind of waves as discussed in Chapter-1. Out of these waves, 
reflected P-wave is the required one i.e. signal and others are unwanted i.e. 
noise in case of common reflection seismology. For better continuity of the 
bed in the seismic section signal to noise ratio should be more i.e. signal 
should have more amplitude compared to noise. Any wavelet contains a 
band of frequencies with different amplitude and phases. Plot of amplitude 
vs. frequency is called Amplitude Spectrum of the wavelet. Signal should 
have wide frequency band i.e. flat amplitude spectrum over wider frequency 
range for better resolution of the bed. Therefore a source should be capable 
of giving us better signal to noise ratio and wide band signal. 


3.1 LAND-SOURCE 

3.1.1 EXPLOSIVE 

The following four parameters decide the wavelet shape and signal to noise 
ratio in the case of explosive energy source. 

(i) Shot hole depth 

(ii) Charge size 

(iii) Geological formation surrounding the charge 

(iv) Shot hole pattern 
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(a) Shot Hole Depth 

Shot hole depth must be greater than the thickness of the low velocity¬ 
weathering layer. At the time of explosion, the force is so large that elastic 
limit of the rock is crossed in the vicinity of the charge, which causes the 
fracturing of the rock material and creation of a cavity. As the material of the 
weathering layer is soft, a cavity of larger radius will be formed causing 
broadening of the pulse width, hence generation of lower peak frequency. 
And since weathering layer is having larger absorption coefficient, most of 
the energy of the wavelet will be absorbed while crossing this layer. Also, 
more nearer the charge to the surface, excitation of earth's upper surface 
would be more, hence generation of more surface waves. Thus shots within 
low velocity unconsolidated layer, cause 

• Generation of surface noise wave. 

• Generation of lower frequency peak. 

• Absorption of useful energy. 

Hence, shot hole depth affects both the frequency and the ratio of signal to 
shot generated noise. Shot depth should be deep enough to give an adequate 
ratio of signal to shot generated noise. 

Experiment for Determination of Weathering Layer Thickness 

For the determination of optimum shot hole depth, first of all, Up-hole and 
Shallow Refraction surveys are done to know the weathering layer 
thickness. 


• Up-Hole Shooting 

Depth and Velocity of weathering and sub-weathering layers can be 
determined by shooting small charge at a number of depths in a shot hole 
and receiving the signals by a detector placed at the top of the hole. Time - 
depth plot is drawn (Fig. 3.1). Depth of different layers can directly be seen 
on depth axis and inverse of the slope of different straight line portion of the 
time - depth curve gives the magnitude of velocities for corresponding 
layers. 
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Fig.3.1: Time-Depth Plot 
• Shallow Refraction Shooting 

At a very short group interval (say five meter) a series of single or bunched 
geophone is planted along a profile and shots are fired from both ends with 
a little explosive at shallow depth say 1 m. Time distance plots are drawn 
for both shots as shown in Figure 3.2. 



Fig. 3.2: time -Distance Plot from both ends (Two Layers Case) 

If earth layers are dipping, we will get two sets of intercepts and slopes one 
each for up dip and down dip case respectively. For calculation purposes 
average of two sets are taken. 

For two layers case 

Sin Ic =Vo /Vi Where, Ic is Critical Angle of incidence 
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Cos Ic = (1-Vo 2 /Vi2)1/2 
Tan Ic=V 0 / (V^-Vo 2 ) 1 / 2 
Therefore, total time 
T=Tab + Tbc + Tcd 
= Z/V 0 Cos Ic + X/Vi + Z/Vo Cos l c 



Fig. 3.3: Ray Path Geometry (Two Layers Case] 

This gives the arrival time for refracted wave as, 

T= X/ Vi + 2Z Cos Ic /Vo 
When X=0, T=Tu 
Therefore, 

Tii = 2Z Cos IC /Vo and 
Z= (Tn/ 2] Vi Vo/V (Vi 2 -Vo 2 ) 

Z can be calculated by substituting the value of Tii. Vo and Vi (from the 
graph) in the above equation. 

In terms of cross-over distance, since at cross-over distance first arrival time 
and refracted arrival times are equal, we get 

Xcross/Vo = Xcross/Vl +2Z{V (Vl 2 - V 0 2 )}/VlV 0 

Z= % V i ( Vl - Vo) / (Vi + V 2 ) U cross 
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Fig. 3.4: Time - distance plot from one end 
(Two layers case] 

For three layers case, ray path geometry is given below 



Fig.3.5: Ray Path Geometry (Three layers case] 
Total time, T = Tab+Tbc+Tcd+Tde+Tdf 
When X=0 , Intercept time 7n is given by 
T| 2 =2Z 0 {V (V 2 2 - Vo 2 )}/V 2 V 0 + 2Zi(V (V 2 2 -Vi 2 )}/V 2 Vi 
Z 1= Vi (T ,2 - 2Z 0 V (V 2 2 - Vo 2 ) / (V 2 V 0 )} (V 2 Vi / V (V 2 2 - V 2 2 )} 

= X/V 2 + 2Z 0 W (V 2 2 - Vo 2 )}/V 2 V 0 + 2Zi{V (V 2 2 - \!^/M 2 Mi 
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Fig. 3.6: Time -distance Plot from one end (three layers case] 

Experiment for determination of Optimum Shot Hole Depth 

Several shots are taken at different shot hole depth below the weathering 
layers keeping charge size fixed. These records are reproduced with 
constant gain, so that all events (signal as well as noise] are clearly visible at 
the level of zone of interest. Peak amplitude of signal & noise is measured for 
a particular channel in all the records at the zone of interest. Peak 
frequencies of the signal are also measured. The depth corresponding to 
record which gives larger signal to noise ratio and larger peak frequency is 
the optimum shot hole depth. 



Fig. 3.7: Amplitude spectrum 
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Amplitude Spectrum of Trace No.l of four records each of depth 22 m, 24 m, 
27 m and 30 m respectively and charge size 2 Kg. is shown in Fig. 3.7. Shot 
hole depth of 22 m shows better performance at higher frequencies in terms 
of bandwidth and peak frequency. 



Fig. 3.8: Traces showing Signal and Noise strengths 

In Fig.3.8 part of monitor record of the same trace number at zone of 
interest from three records of different depths are shown. Trace from record 
of 27 m depth shows the highest signal to noise ratio. 



Shot Hole Depth 9m 
Fig. 3.9 


Shot Hole Depth 20m 
Fig. 3.10 


Shot Hole Depth 27m 
Fig. 3.11 
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Fig. 3.9, Fig.3.10 and Fig. 3.11 are the field records for 9m, 20m and 27m 
shot hole depth respectively. 

(B) Charge Size 

As discussed in the previous section that blasting of explosive creates a 
cavity in the surrounding medium. Radius of the cavity determines the pulse 
width of the wavelet. Since the smaller charge will create a cavity of shorter 
radius, thus produces shorter pulse width and higher frequency. It has been 
observed that the time duration and amplitude of the pulse produced by 
explosion is proportional to M 1 / 3 , where M is the mass of the spherical 
charge. Smaller charge will give larger frequency. 

Charge size should be determined by the ratio of signal to ambient noise. 
The acceptable ratio at the target depth varies from 1:1, when maximum 
resolution is required provided multiplicity of the recorded data is very high 
and to 10:1 when high quality data is required in the case of low multiplicity. 

Experiment for determination of Optimum Charge Size 

To determine the optimum charge, shots with different charge size but all 
having optimum depths are taken and signal to noise ratio at the level of 
zone of interest is computed. Peak frequency of the signal is also computed. 
The record with lowest charge size which gives acceptable signal to noise 
ratio corresponds to optimum charge size. 

Formation Surrounding the Charge 

Formation surrounding the charge also plays a major roll to decide the 
required wavelet shape. Charge blasted in medium of sand and gravel gives 
low frequency peak and lower amplitude of the signal, while the medium 
like hard clay gives higher frequency peak as well as higher amplitude of the 
signal. Therefore this factor should also be considered while deciding the 
optimum depth. 

(C) Shot Hole Pattern 

In the area, where surface noise cannot be suppressed by other means, 
pattern of shot holes just like geophone pattern may be used to reduce the 
effect of ground roll. If apparent wavelength of the prominent noise 
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traveling along the profile is ‘ A,n ' then shot hole spacing'd ' of n holes 
pattern will be 

d = A, N / n 

If in any area, off-line noise is a problem, then areal pattern of shot holes 
may give satisfactory results. If there are two prominent surface noise 
having different wavelengths, then one may be suppressed with receiver 
array and other with the help of shot array. 

3.1.2 AIRGUN 

Airgun is a surface energy source but just like explosive it is also an 
impulsive source. With the help of 'micro processor based' firing control 
system any number of airgun can be fired simultaneously with a standard 
deviation of the order of 200 microseconds. Since it is low energy surface 
source, number of air guns in use, number of pops for each air gun and air 
gun pattern are the three important parameters to decide signal to noise 
ratio. The frequency spectrum of the land air gun is dependent on a number 
of factors, such as chamber volume, air pressure, surface area of the base 
plate and surface and near surface condition of the earth. A typical source 
signature in time and frequency domain produced by airgun is shown in 
Fig. 3.12 & 3.13 respectively. 



Fig. 3.12: Typical Source signature in time domain 
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Fig. 3.13: Typical Source signature in frequency domain 


(A)Number of Air guns & Number of Pops 

Field experiment have shown that the law of linear superposition holds very 
closely in the case of multiple air gun operating simultaneously i.e. all 
synchronized with control signal. Therefore any number of air gun can be 
used simultaneously. Number of air guns and number of pops depend on the 
signal to random noise ratio, we want to achieve. Since air gun is a surface 
energy source it produces very low energy reflection as compared to surface 
noise. Therefore, more number of air gun or more number of pops per air 
gun to yield a single record will produce signal of better energy. If there are 
'n' air guns in use then signal to random noise ratio will be 'n' times more 
than the single air gun. If'p' number of pops are used per air gun, then signal 
to random noise ratio after vertical stack will be n.p 1/z times more than the 
single air gun with single pop. For, example, if desired maximum signal to 
noise ratio be 18 dB (dB=201og n.p 1/z ), the relative number of pops Vs. 
number of air gun would be 

1 gun 64 pops 

2 gun 16 pops 

3 gun 7 pops 

4 gun 4 pops 

Figure 3.14 shows the relationship between signal to noise ratio in dB and 
number of pops for different number of guns. 
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Fig. 3.14: S/N ratio in dB and number of pops for different number of guns 

(B) Air-Guns Pattern 

The increase in number of pops or number of air guns only improves the 
signal to random noise ratio. To reduce ground roll or any coherent noise 
suitable air gun pattern may be designed. Just as receivers array, a linear or 
tapered array of air guns may be used, depending upon the noise 
wavelength and degree of attenuation we want. Criterion of design of source 
array parameter is the same as receiver array parameters. Only limitations 
is the number of air guns available and number of pops we can afford to use. 
If is the prominent noise wavelength and 'n' is the number of source 
elements we want to use, then element spacing would be 

d = A, n / n 

Let us suppose that we have four air guns A, B, C, D and we want to use 12 


elements source pattern with spacing'd'. 

1 2 3 

4 

5 

6 7 8 9 10 11 12 

A 

B 


C D 

A 


B 

C D 

A 



BCD 


In this case airgun A will be kept at element position 1, airgun-B at position 
4, airgun-C at position 7 and airgun-D at position 10. After taking one shot at 
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this position, A will move to position 2, B to 5, C to 8 & D to 11 and second 
shot will be taken. After taking second shot A will move to 3, B to 6, C to 9 
and D to 12 and third shot will be taken. When all the 3 shots are vertically 
stacked, it will produce the effect of 12 elements source array. 


1 2 

A B 

A 


3 

C 

B 

A 


4 

D 

C 

B 


5 

D 

C 


6 

D 


If we want to use tapered source array 1-2-3-3-2-1 then, the entire four air- 
guns are kept at positions 1 to 4 in sequence and shot is taken. Then each 
airgun moves one position ahead i.e. A goes to position 2, B to 3, C to 4 and D 
to 5 and second shot is taken. Lastly A goes from 2 to 3, B from 3 to 4, C from 
4 to 5 and D from 5 to 6 and third shot is taken. Vertical stack of these three 
records will give the result of the array 1-2-3-3-2-1. In this way any desired 
array can be designed to attenuate ground roll, as done in the case of 
receiver arrays. 

(C) Synchronization Test of Air gun 

All the air guns in use must fire simultaneously within permissible limits. 
Otherwise there will be a phase shift between the signals produced by each 
air gun and combined effect will be loss of high frequency and amplitude as 
well. Therefore, this test must be performed regularly. 

3.1.3 VIBRATOR 

In the Vibroseis (trade name) system a signal of long duration and changing 
frequency (called sweep) is applied to the earth. The reflections received are 
delayed replicas of the input sweep signal, but as the sweep is long, 
reflections overlap one another unlike the reflections from explosive and 
airgun sources. This prevents determination of arrival time of each 
reflection. These overlapping reflections must be separated by compressing 
them to their impulse equivalents. This is accomplished by cross correlation 
of reflected signal with the pilot sweep signal. The pilot signal will locate the 
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replica of it and produce a correlation wavelet at the place of reflection. Even 
though these correlation wavelets are not the actual signals that traveled 
into the earth, they however have same spectral properties as the one that 
has traveled. Theoretical auto- correlation wavelet of sweep is called 
KLAUDER WAVELET. Since the actual cross-correlated wavelet resembles 
the Klauder wavelet, the study of the characteristics of Klauder wavelet is 
very helpful for deciding the input sweep. 

(A) Characteristics of Klauder Wavelet 

The idealized wavelet is the auto-correlation function of the linear swept 
frequency signal input to the ground. When the duration of the input signal 
is greater than one second, the autocorrelation is closely approximated by 
the expression 

K(t]=T [{ Sin n A (1- t/T) t} / IT A t ]Cos2ITfot (3.1) 

Where, 

T= duration of the sweep in seconds, called sweep length. 

A = f 2 - fi = bandwidth 
f 2 =limiting high frequency of the sweep, 
fi = limiting low frequency of the sweep. 
fo = mid-frequency of the sweep =ffi+f 2 )/2 
K(t) = Klauder wavelet. 

The shape of the auto-correlation wavelet depends on the mid frequency fo 
and the ratio of the terminal frequency (f 2 /fi) of the input sweep. 

The shape or character of the side lobes depends on the ratio of terminal 
frequencies. The ratio of the amplitude of the central peak to the amplitude 
of the first side lobe is a function of the ratio of the highest to the lowest 
frequency of the input sweep. The zero crossing interval for the central peak 
of the autocorrelation wavelet is half of the period 'To' of the mid frequency 
‘fo’. The Klauder wavelet has a low frequency underlying the high frequency 
jitter of the side lobes. This low frequency is approximately the same as the 
lowest frequency of the input sweep and the superimposed high frequency 
jitter is the same as the highest frequency end of the input signal. Variation 
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of autocorrelation shape with respect to frequency ratio (f 2 /fi), mid 
frequency (fo] and band- width (f 2 -fi) of input signal is shown in Fig, 3.15. 

On the left are two octave signals. The wavelets have essentially identical 
shapes except for the time scale. 



The higher frequency signal at the bottom is compressed in time. But the 
shapes i.e. the relative amplitude of each successive up or down excursions 
are the same. 

In the center all the wavelet are having the same central frequency 20 Hz. It 
is quite clear that the smaller the ratio of high to low frequency the more 
legs there are in the wavelet and the Klauder wavelet appears to ring more. 

On the right, all the Klauder wavelet in Fig. 3.15, have same bandwidth (f 2 -fi) 
24 Hz. At the top is the wavelet for 4-28 Hz signal, which has a bandwidth of 
about 2.8 octave and has a good wavelet shape. Below that is 8-32 Hz i.e. 2 
octave signal, which also has a good wavelet shape. The 16-40 Hz. 1.3 octave 
wavelets has significant amplitude on the first side lobe as compared to 
central peak. When the frequency ratio drops, side lobes amplitude increase 
and wavelet shows more number of legs. 
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(B) Optimization of Vibroseis Input signal 
Selection of best Klauder Wavelet 

Our basic objective is to choose an input signal that would produce a 
seismic wavelet (Klauder wavelet] which would give the best resolution of 
the geologic problem. To find this wavelet, we make a spike synthetic 
seismic trace that represents the geology of the interested horizons. Then 
we convolve the spike synthetic trace with Klauder wavelet having different 
input frequency ranges and different centre frequencies to find the wavelet 
that give the best resolution of the geologic objective. 

The linear-sweep frequency signal having autocorrelation function same as 
the best Klauder wavelet is, theoretically, the best signal to use in the field to 
fulfill the sub-surface geologic objective. If there are a number of different 
geologic objectives, each requiring a different ideal wavelet shape, then we 
must compromise in specifying the input signal frequencies. Generally, we 
would concentrate on the geologic zone, with the greatest economic value. 

Selection of Optimum Low Frequency end of Sweep 

Having determined the theoretically optimum input signal, the next step is 
to make field tests to determine the practically optimum input signal. 
Supposing the 3-octave input signal 10-80 Hz gives the best synthetic 
correlogram, 10-80 Hz input, was decided on the basis of synthetics. Several 
alternate signals, whose end frequencies are near about 80 Hz and beginning 
frequencies are different such as wavelets 12-81 Hz, 14-79 Hz, 16-76 Hz, 17- 
81 Hz, 20-80 Hz, 22-74 Hz, 24-80 Hz, 28-80 Hz are used and separate record 
with full spread is taken for these different input signal. 

Suppose the dominant frequency of strong noise wave in the area is 14 Hz, 
then an input signal with ‘fT greater than 14 Hz should produce a better 
signal to noise ratio. If noise analysis from several test sites showed that the 
surface wave frequencies ranged between 10 and 16 Hz, then 17-81 Hz 
input signal would be a likely compromise between the theoretically best 
10-80 Hz signal for defining the desired geologic horizons and the require¬ 
ment of minimizing surface wave interference. 
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Selection of Optimum High Frequency end of Sweep 

Suppose the record shows a strong wind noise in 70-150 Hz band. Then 81- 
150 Hz interference will be reduced by cross-correlation because it is 
outside the previously decided, 17-81 Hz input signal band. But 70-81 Hz 
portion could be so strong that it would mask the high frequency portion of 
the reflection. To choose the optimum value of the high frequency end of the 
sweep, different records with different high frequencies end and having 
fixed low frequency end are taken with full spreads. Let us suppose that 
available sweeps are 17-68, 17-48, 17-40 17-34 Hz which covers the 
bandwidth ratio range from 2.25 to 1 octave. Then good choice may be 17- 
57 Hz or 17-68 Hz input signal for attenuation of high frequency noise. 

The above field tests should be repeated at several locations within the 
prospect, especially if there are appreciable differences in near surface 
geology. 

Selection of Optimum Sweep Length 

Sweep length, of the signal is important for two reasons. Firstly, the 
amplitude of correlation wavelet is larger for larger sweep length as is 
evident from equation (3.1). Second is to keep the correlated harmonics of 
low frequency, outside the zone of interest. The second point can be 
understood in the following way. Ideally the input to the ground is a copy of 
the pressure applied to the base plate but in fact crushing and compaction of 
the surface material result in the input to the ground varying nonlinearly 
with the pressure exerted by the vibrator, this introduces harmonics. This 
effect is less serious for down sweeps but up-sweeps are much easier on the 
equipment. So generally up-sweeps are preferred. Let us suppose that input 
to the ground is 20-80 Hz signal with sweep length 6 sec. Therefore, rate of 
change of frequency is 10 Hz per sec. and signal after 2 sec. will attain the 
frequency of 40 Hz. First harmonics of 20 Hz is also 40 Hz, so there is chance 
of getting harmonics of 20 Hz correlated with the proper fundamental at 2 
sec. and producing a ghost of higher frequency. If this is the case, we must 
keep this ghost outside the zone of interest. This can be done by increasing 
the sweep length. Suppose that the sweep length is increased to 12 sec. then 
rate of change of frequency in the input sweep would be 5 Hz/sec it would 
take 4 sec for attaining the frequency of 40 Hz input frequency. If our main 
zone of interest is 3 sec, then correlated harmonics i.e. ghost will fall outside 
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the zone of interest. In this way we can get rid of the ghost from the zone of 
interest. Effect of different sweep lengths is shown in Fig. 3.16. 

Sweep Sweep Sweep Sweep Sweep 

Length: 7 Sec. Length: 6 Sec Length: 5 Sec. Length: 4 Sec Lengths: 3 Sec 



Fig.3.16: Effect of Sweep length on Ghost Harmonic 


In this figure event GG1 is the ghost of event PP1. As the sweep length 
increases event GG1 shifts towards higher times of the record and ultimately 
goes outside the zone of interest. 

Selection of Optimum number of Sweeps 

Since Vibroseis is a low energy source, several sweeps of input signals are 
required to get better signal to noise ratio after vertical stack of all the 
sweeps. Signal to random noise improvement is proportional to square root 
of the number of sweeps. Records with different number of sweeps are 
taken and the record which gives acceptable level of signal to noise ratio 
may be considered as optimum from the point of view of number of sweeps. 

Selection of Vary Sweeps 

Some times limiting frequencies fl and f2 are varied for each sweep to fill 
the notches in the amplitude spectrum of the wavelet and to reduce the side- 
lobes in correlation wavelet. For example instead of one fixed sweep 14-80 
Hz it may be better to have six sweeps of frequency bands 14—56Hz, 16- 
64Hz, 17—63Hz. 18—72Hz, 17—96 Hz, 20—80Hz, which contains less 
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number of side lobes than the individual sweep. To decide optimum vary 
sweep, records with different combination of sweeps are taken. Shape 
and amplitude spectrum of the cross-correlated wavelets corresponding to 
each combination are compared. The combination of sweep, which gives less 
number of notches in the amplitude spectrum and less number of side lobes 
is the optimum vary sweep. 

Selection of Vibrator Pattern 

Law of linear superposition holds good in the case of multiple vibrators 
operating simultaneously. Therefore any number of vibrators can be used at 
a time, all synchronized with control signal. Since length of vibrator input 
signal is very large (typically 5 to 12 seconds long] as compared to impulsive 
source, compressional and shear head waves and Rayleigh type surface 
waves completely cover almost all of the reflections. Therefore, if we record 
on digital tape with a broad frequency band at full scale, it is the amplitude 
of shear head wave and ground roll that determine the gain in the amplifier 
rather than the amplitude of the reflections that we would like to record. 
Therefore, the vibrators should be spaced to provide a source pattern that 
would attenuate wavelengths in the head waves so that the seismometers 
would receive a lower amplitude head wave and thus improving the signal 
to noise ratio, This reduction in amplitude of the head waves is one of the 
most important techniques for achieving maximum signal to noise ratio in 
vibrator work. 

The head waves will ordinarily have one velocity, the sub-weathering 
velocity unless there is enough distance between the source and the farthest 
receivers to allow several refraction branches of multiple sub-weathering 
layers to be recorded. This velocity is divided by the dominant period of 
head wave to give the apparent wavelength of head-wave. A suitable pattern 
of vibroseis may be designed to attenuate this head wave as we decide for 
the receiver pattern. 


3.2 MARINE SOURCE 

Choice of energy source is a major part of planning for proper data 
acquisition. The size and complexity of the source have direct bearing on the 
cost of data acquisition. In general, larger sources produce more usable 
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seismic energy than smaller source, hence may improve S/N of the data. On 
the other hand larger source (say 100 bar-m] may be more expensive than 
the smaller source (say 25 bar-m] this proper source specification is 
essential for proper data acquisition with minimum cost. 

The vertically traveling far-field signature of air gun array produces a 
quantitative measure of array performance. A typical gun signature has been 
shown in Fig 3.17. 

The large positive peak in the signature is called the direct or primary 
arrival. It corresponds to the increase in pressure that occurs in the water as 
the compressed air first explodes out of the guns, the primary arrival is 
followed by a large negative peak, the surface reflected or ghost arrival. The 
time delay between the primary and the ghost arrivals depends on the 
source depth and hydrophone position. The ghost arrival is a negative peak 
because the reflection coefficient of the water air interface is negative. After 
the above two main arrivals there are a sequence of a small variation in 
pressure called the residual bubble pulse. These features correspond to the 
pressure variation that occurs in the water due to expansion/compression 
of the air bubbles. One of the main goals of air gun array design is to make 
the bubble pulses as small as possible compared to the primary arrival. The 



Fig 3.17: Typical gun signature in time domain 








86 


N 

X 


I 

M 

CL 


it 

> 

■a 

_2 

a 

u 

cc 

T3 



Fig 3.18: Typical gun signature in frequency domain 

spectrum has a large scale variations consisting of peaks and notches. 
These correspond to frequency at which there is constructive and 
destructive interferences respectively between the primary and ghost 
arrival and the smaller scale variation in the spectrum between 10-100 Hz 
is caused by the bubble oscillation. 

3.2.1 SPECIFICATION OF SOURCE SIGNATURE 

There are two specifications which are commonly used, (i) giving the 
signature strength and (ii) describing the signature shape; the two aspects 
can be formulated either in time domain or frequency domain. 


(A) SIGNATURE STRENGHT IN TIME DOMAIN 

A far field hydrophone detects pressure wave during the signature test but 
because of spherical spreading the magnitude of the pressure wave depends 
on how far the hydrophone is from the array. Therefore the source strength 
is quoted in terms of pressure-meter [bar-m] unit for example 25 bar-m 
means that if an array were a point source and a hydrophone were one 
meter away then the hydrophone would detect a pressure of 25 bar, if the 
phone were 10 m away it would detect a pressure of 2.5 bar (spherical 
spreading correction l/r].The peak strength is simply the amplitude of the 
primary arrival. The proper SI pressure unit used by others in underwater 
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acoustics is the micro-Pascal, equivalent to 10-11 bar. The peak to peak 
strength (PTP] also include the amplitude of ghost arrival and thus is the 
distance between the two major extremes of the signature in the Fig 3.17 
the peak to peak strength is given by the separation between dashed lines A 
and B. 

(B) SIGNATURE SHAPE IN TIME DOMAIN 

Signature shape is represented in time domain by the primary to bubble 
ratio (PBR], This is computed by dividing the PTP signature strength by the 
peak to peak strength of bubble pulse. In the Fig 3.17 

PBR = (separation of line A and B] / (separation of line C and D] 

Air gun array with high PBR is called tuned air gun arrays. 

(C) SIGNATURE STRENGHT IN FREQUENCY DOMAIN 

In frequency domain, signature strength is specified by absolutely 
calibrating the signature amplitude spectrum in dB relative to IpPa / Hz at 
lm, as shown in Fig 3.18 it is impossible to compare performance of two 
sources in frequency domain if source spectra are referenced to 0 dB. 

(D) SIGNATUR SHAPE IN FREQUENCY DOMAIN 

In frequency domain, signature shape is specified by describing the 
smoothness of the signature's amplitude. For example in the Figure 3.18, we 
could say that between 10 and 100 Hz the spectrum is flat within ± X dB 

(E) EFFECT OF BANDWIDH ON TIME DOMAIN SPECIFICATIONS 

Let us suppose the signature has 1msec sampling interval and consequently, 
a Nyquist frequency of 500 Hz and as a practice 375 Hz anti-alias filter was 
applied. Therefore signature has been band-limited with 375 Hz high cut. It 
is observed that as the high cut frequency decreases, the PTP and PBR of an 
array also decreases. Time domain specifications are also affected by 
application of low cut, partly due to change in the low frequency content of 
the signature and partly due to the phase change inherent in the instrument 
filter. Thus the signature from different arrays cannot be easily compared by 
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time domain specifications if they have either different low-cut or high-cut 
filters. 

(F) EFFECT OF ARRAY PARAMETERS ON SPECIFICATIONS 


PARAMETERS 

PTPSTRENGTH 

PBRSTRENGTH 

Number of guns 

Proportional to number 
of guns 

No direct relation 

Gun volume 

Proportional to the cube 
root of gun volume 

Increases as volume 

increases 

Gun depth 

Increases at first, then 
decreases as depth 
increases 

Decreases as gun depth 
increases 

Initial firing 
pressure 

Increases as firing 
pressure increases 

Increases as firing 
pressure increases 

Port area 

Proportional to square 
root of port area 

Increases as port area 
increases 

Port-closure 

pressure 

No effect if closure 
pressure < 0.5 of firing 
pressure 

Decreases as port-closure 
pressure increases 

Gun 

synchronization 

Decreases as spread in 
firing times increases 

Decreases as spread in 
firing times increases 

Horizontal gun 
separation 

Decreases as guns move 
closer together 

Increases dramatically 
when gun bubbles 
coalesce 
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CHAPTER-4 


RECEIVER AND SOURCE ARRAYS 


4.0 INTRODUCTION 

W hen earth is excited either with surface sources or explosive sources 
blasted in a shot hole, surface waves are also generated in addition to 
the P-wave. This surface wave travels with velocity around 300 - 600 m /sec. 
and generally has frequency around 5 Hz to 20 Hz. Amplitude of surface 
wave is normally very large compared to the amplitude of reflected P-wave. 
Therefore, it is necessary to suppress this wave before entering into the 
recording instrument. Since surface wave has low frequency and low 
velocity compared with the signal, the apparent wave number of surface 
wave is smaller than that of reflected P-wave. Therefore suitable spatial 
filter is designed to discriminate the waves according to wave numbers. This 
can be done in the field by designing proper arrays of sources / receivers, 
instead of using a single source/receiver. Design of array depends on the 
characteristics of surface waves and requirement of resolution in a 
particular area. 


4.1 NOISE EXPERIMENT 

To find wavelengths and frequencies of different coherent waves in the area, 
special experiment is carried out in the field. Full spread is laid on an even 
surface keeping single geophone string bunched at very small group interval 
(generally 5 m], Length of the spread will be 5 [n-1] m, where n is number of 
channels and 5 m is group interval as assumed. Therefore, for 96 channels 
spread length will be 475 m. First shot is taken with optimum source 
parameters as already decided at an offset of one group interval (i.e. 5 m] 
from the nearest group. Therefore, for the first shot, nearest channel (1 st 
channel] is at 5 m and farthest channel (last channel] is at 480 m. Then for 
the second shot keeping shot location fixed at the same position, spread is 
moved by a spread length (i.e. 475 m] so that now first channel is at 480 m. 
and last channel is at 955 m. Similarly 3 rd , 4 th and 5 th shots are taken at the 
same shot location and spread is moved each time by one spread length as 
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shown in Fig. 4.1. Almost the same result can also be obtained by keeping 
the spread fixed at the same location and moving the shot point locations 
each time by one-spread lengths, (i.e. 475 m in this case] as shown in Fig. 
4.2. Both these types are called in line wave-test. Sometimes cross-line 
wave-test is also performed to receive near surface scattered wave coming 
from off-line direction. Spread in cross-line is kept perpendicular to the in¬ 
line spread. 


FIXED SP. First Spread 

Third Spread 

SP. CH-1 CH-96 

CH-1 

CH-96 

A # # 

£ 

# # 

CH-1 

CH-96 

CH-1 CH-96 

Second Spread 

Fourth spread 


Fig. 4.1: Noise Test (Moving Spread] 


SP-4 

SP-3 

SP-2 

SP-1 

FIXED SPREAD 

CH-1 C H-96 

A 

A 

A 

A 
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Fig. 4.2: Noise Test (Moving Shot Points] 



Fig. 4.3: Field Record for Noise Test (Moving Shot Points] 
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Different record shows the behavior of surface waves for different distances 
from the shot point. First record shows the behavior of surface wave within 
5-480 m. Second record shows the behavior of surface wave within 480 - 
955 m. Third record show the behavior of surface wave within 955-1430 m. 
Fourth record shows the behavior of surface wave within 1430-1905 m. 
Fifth record shows the behavior of surface wave within 1905 - 2380 m. Sixth 
record shows the behavior of surface wave within 2380 - 2855 m. In this 
way behavior of surface wave is recorded up to far offset distance, which has 
been decided to use in regular production work. A field record for noise test 
is shown in Fig.4.3. 


4.2 ANALYSIS OF NOISE RECORDS 


Time-distance (t-X] graph is plotted for each wave observed in the field 
record. Each wave will align along different straight line; slope of each 
straight line gives the apparent velocity of the corresponding wave as shown 
in Fig.4.4. Time period of the wave can be found by measuring the time 
difference between two consecutive peaks or troughs where the wave shape 
is clear and with least interference. Inverse of the period gives the 
frequencies of the corresponding waves. In the same way frequencies and 
apparent velocities of signal at different offsets and times can also be 
measured from the record. Amplitudes of signals and noises can also be 
measured at different offsets & times. 



Fig.4.4: Time-distance (t-x) plot for noise 
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NOISE-A: f=14 Hz, V=440 m/s 

NOISE-B: f= 9 Hz, V=370 m/s 

NOISE-C: f= 5 Hz, V=300 m/s 


4.3 ARRAY THEORY 

Attenuation of noise especially surface waves is the most important factor 
for enhancement of signal to noise ratio S/N in the field. Coherent noise 
(traveling events which can be followed from trace to trace for substantial 
distances along the profile] can be suppressed by applying principle of 
directivity similar to those used in designing antenna for receiving radio- 
waves. Detectors are grouped in arrays for two purposes, first to attenuate 
random noise or incoherent noise (scattered wave from near surface 
irregularities] arriving from any direction and second to attenuate coherent 
noise arriving from specified direction. Cancellation effect of random noise 
is proportional to VN where, N is the number of elements in the array. Thus 
to obtain maximum suppression of incoherent noise one would therefore 
use as many elements per trace as possible. 

4.3.1 RESPONSE OF UNIFORM LINEAR ARRAY 

(Equal Weight - Equal Spacing] 

For linear array of equal weights with number of elements N (sources or 
receivers] at spacing ' d ’ and each element having sensitivity equal to ao,the 
resultant output of the array due to the action of a harmonic wave of angular 
frequency eo will be 

Ao = ao exp (icot] + ao exp (icot-i®] + ao exp (icot -2i®] +. 

+ ao exp [ieot-i(N-l]®] 

Where ® = (2 lid Sin 0 / A.] is the phase difference between two elements 
spaced at a distance'd ' for the wavelength X and 0 is the angle which 
traveling wave makes with the normal to array direction. Therefore, 

Ao = ao [1 + exp (- i®] + exp (-2 i®] +... + exp {- i (N -1]} ®] exp (i co t] 

Ao = ao [{1- exp (-i N®]} / (1- exp (- i®]]] exp (i co t] 

Therefore, amplitude of the output wave will be given by 
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Ao = ao [{1- exp (- i N ®]}{1- exp (i N ®]} / 

{1- exp (- i ®]} {1- exp (i O]}] 1 / 2 

Ao = ao [{2- exp (- i N®] - exp (i N ®]} / {2- exp(- i ®] - exp (i ®]} ] y2 

Ao = a 0 [{1- Cos N® } / {1- Cos ® }]V2 

Ao = a 0 [{Sin (N® / 2]} / { Sin (® / 2)}] 

Ao = ao [{Sin (N tc d Sin 0 / X~)}/ {Sin { n d Sin 0 / X ]}] 

If the wave is traveling along the direction of the array, then 0 = n /2 and 
therefore, 

Ao = ao [{Sin (N n d /X)} / {Sin ( n d / X )}] [4.1] 

If N elements are bunched at one place, the resultant amplitudes will 
be B 0 =N a 0 

Therefore relative amplitude response of the array with respect to the 
bunched elements will be 

Ar = Ao / Bo 

Ar = [{Sin { N 7i d Sin 0 / X ]} / N{Sin { n d Sin 0 / X ]}] 

And if 0 = tc / 2, 

Ar = [{Sin (N n d / / N{Sin { tc d / X ]}] 

Since the power is proportional to the square of amplitude, the Relative 
power response in terms of decibel will be given by, 

R= 10 log [Ar] 2 

R = 20 log [{Sin (N n d Sin 0 / X ]} / N{Sin { n d Sin 0 / X )}] dB 
And if 0 = n / 2, 

R = 20 log [{Sin (N nd / ^ ] } / N{Sin [itd/l]}/N] dB (4.2i) 

Ao/ao for different number of elements has been calculated for different 
values of A/X and tabulated in Table 4.1. With the help of this table array 
response curve in dB for different number of elements can be drawn by 
dividing the values of Ao /ao by number of elements N, taking logarithm and 
multiplying by twenty, which gives the value of R in dB. 
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Fig.4.5: Response curve for 6-element uniform array 

Response curve for a uniform linear array of six elements is shown in the 

Fig. 4.5., it is observed that, 

(a] The response curve is symmetrical about d/^= 0.5 which corresponds 
to the nyquist spatial wave number. 

(b] The notches occur at wave number 1/Nd, 2/Nd, (N-1]/Nd. 

Where N=2, 3, 4,. 

(c] Secondary peaks occurs at wave number 3/2Nd, 5/2Nd, (2N—3]/ 2Nd, 

Where N = 3, 4, 5, 6. 

(d] First alias peak occurs at wave number 1/d. 

(e] Relative amplitude at the bottom of the envelope 'E / is proportional to 
1/N 

(f] Large number of detectors should be used for optimum attenuation at 
the center of the lobe 
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Table 4.1 


Amplitude Response Ao = ao [{Sin (N7id / A.)} / (Sin ( n d / 1 )}] (for linear array) 


d/1 

NO . OF ELEMENTS 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

0.01 

1.999 

2.996 

3.99 

4.98 

5.966 

6 . 9fl5 

7.917 

8.882 

9.838 

10.78 

11.72 

0.02 

1.996 

2.984 

3.961 

4.921 

5.863 

6.781 

7.672 

8.534 

9.361 

10.15 

10.902 

0.04 

1.984 

2.937 

3.844 

4.69 

5.462 

6.148 

6.737 

7.219 

7.588 

8.837 

7.963 

0.06 

1.965 

2.860 

3.653 

4.317 

4.829 

5.169 

5.326 

5.295 

5.076 

4.677 

4.112 

0.08 

1.937 

2.753 

3.395 

3.824 

4.013 

3.95 

3.638 

3.098 

2.364 

1.48 

0.504 

0.1 

1.902 

2.618 

3.078 

3.236 

3.078 

2.618 

1.902 

1 

0 

1 

1.902 

0.12 

1.860 

2.458 

2.711 

2.584 

2.093 

1.309 

0.34 

- 0.676 

- 1.597 

- 2.294 

- 2.688 

0.14 

1.810 

2.275 

2.307 

1.9 

1.131 

0.147 

- 0.865 

- 1.712 

2.234 

2.33 

- 1.983 

0.16 

1.753 

2.072 

1.878 

1.22 

0.26 

- 0.764 

- 1.600 

- 2.039 

- 1.974 

1.421 

- 0.516 

0.18 

1.689 

1.852 

1.438 

0.577 

- 0.464 

- 1.360 

- 1.833 

- 1.735 

- 1.097 

0.117 

0.899 

0.2 

1.668 

1.618 

1 

0 

- 1.0 

- 1.618 

- 1.618 

- 1.0 

0 

1 

1.618 

0.22 

1.541 

1.375 

0.578 

- 0.485 

- 1.325 

- 1.556 

- 1.074 

- 0.099 

0.922 

1.52 

1.42 

0.24 

1.458 

1.126 

0.138 

- 0.859 

- 1.435 

- 1.233 

- 0.363 

0.704 

1.389 

1.322 

0.538 

0.26 

1.369 

0.874 

- 0.172 

- 1.110 

- 1.348 

- 0.735 

0.341 

1.202 

1.305 

0.584 

- 0.505 

0.28 

1.275 

0.625 

- 0.478 

- 1.234 

- 1.096 

- 0.163 

0.888 

1.295 

0.763 

0.323 

- 1.174 

0.3 

1.176 

0.382 

- 0.727 

- 1.236 

- 0.727 

0.382 

1.176 

1 

0 

- 1.0 

- 1.176 

0.32 

1.072 

0.148 

- 0.913 

- 1.126 

- 0.295 

0.811 

1.163 

0.436 

- 0.696 

- 1.182 

- 0.571 

0.34 

0.964 

- 0.072 

- 1.033 

- 0.923 

0.143 

1.061 

0.879 

- 0.214 

- 1.085 

- 0.832 

0.284 

0.36 

0.852 

- 0.275 

- 1.086 

- 0.650 

0.532 

1.103 

0.407 

- 0.757 

- 1.051 

- 0.139 

0.933 

0.38 

0.736 

- 0.458 

- 1.073 

- 0.332 

0.829 

0.942 

- 0.135 

- 1.042 

- 0.632 

0.576 

1.056 

0.40 

0.618 

- 0.618 

- 1.0 

0 

1 

0.618 

- 0.618 

- 1.0 

0 

1 

0.618 

0.42 

0.497 

- 0.753 

- 0.872 

0.319 

1.03 

0.193 

- 0.934 

- 0.658 

0.607 

0.96 

- 0.129 

0.44 

0.375 

- 0.860 

- 0.697 

0.598 

0.921 

- 0.253 

- 1.007 

- 0.128 

0.968 

0.49 

- 0.784 

0.46 

0.251 

- 0.937 

- 0.486 

0.815 

0.69 

- 0.642 

- 0.851 

0.429 

0.959 

- 0.189 

- 1.006 

0.48 

0.126 

- 0.984 

- 0.249 

0.953 

0.369 

- 0.907 

- 0.686 

0.846 

0.589 

- 0.772 

- 0.686 

0.50 

0 

- 1.0 

0 

1 

0 

- 1.0 

0 

1 

0 

-1 

0 
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Table 4.2 


No of 

elements 

1st 

Peak 

2nd 

Peak 

3rd 

Peak 

4th 

Peak 

5th 

Peak 

6th 

Peak 

7th 

Peak 

8th 

Peak 

9th 

Peak 

10th 

Peak 

3 

-9.4 










4 

-11.4 

-11.4 









5 

-12.1 

-14.0 

-12.1 








6 

-12.6 

-15.3 

-15.3 

-12.6 







7 

-12.8 

-16.0 

-16.9 

-16.0 

-12.8 






8 

-13.0 

-16.5 

-17.9 

-17.9 

-16.9 

-13.0 





9 

-13.1 

-16.8 

-18.5 

-19.1 

-18.5 

-16.8 

-13.1 




10 

-13.14 

-17.0 

-19.0 

-20.0 

-20.0 

-19.0 

-17.0 

-13.24 



11 

-13.18 

-17.15 

-19.33 

-20.47 

-20.83 

-20.47 

-19.32 

-17.15 

-13.18 


12 

-13.24 

-17.27 

-19.57 

-20.90 

-21.51 

-21.51 

-20.9 

-19.57 

-17.27 

-13.24 

13 

-13.25 










14 

-13.3 










20 

-13.38 










50 

-13.45 










100 

-13.46 











Table 4.2 shows the values of response 'R' at secondary peaks of response 
curves in dB for uniform arrays of different elements, indicating that, as the 
number of elements increases the values of R at first and successive peaks 
decreases. But we see there is a lower limit to this value at first peak, which 
is approximately - 13.465 dB. 

The array output lags behind that of the first geophone by an angle (Fig. 4.6] 
a = (N-l}<|) / 2 
Where, § = 2 n d Sin 0/ X, 
d = Element spacing, 

N = Number of elements. 
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The width of the principal lobe for uniform array depends on the value (N.d], 
which is one geophone spacing greater than the actual length of array; (N.d] 
is called effective array length. The region between the points where the 
response is down by 3 dB is called the reject region, because power in this 
region is less than half of the peak power. 

Generally longest noise wavelength & shortest signal wavelength 
component, which is the shallowest reflection of interest on the longest 
offset trace (highest frequency at lowest apparent velocity] overlap 
considerably. Array should be designed in such a way that it should not 
affect the signal adversely. By measuring the apparent velocity and 
estimating the maximum frequency component we can estimate the 
wavelength of this reflection component. The array, which has been 
designed to attenuate longest wavelength noise, must not attenuate the 
signal of shortest wavelength of interest by more than 3 dB, if maximum 
possible resolution is desired. The noise wavelength longer than this may be 
eliminated by suitable choice of low cut filter of recording instrument, if 
necessary. 

4.3.2 RESPONSE OF ELEMENTS WEIGHTED TAPERED ARRAY 

(Unequal Weights & Equal Spacing] 

Array where different elements have different weights are called element- 
weighted tapered array. The main lobes of response curve are widened but 
the attenuation in the reject region is generally large compared with a linear 
array with the same overall array length. For the areas where S/N ratio is 
good, a linear uniform array pattern is sufficient. Weighted arrays are used 
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in the areas where high amplitude noise affects the S/N ratio much, and 
more attenuation of noise is required. 


(A) Method Of Calculation Of Response Curve 

For the sake of calculation any weighted array may be considered to be the 
combination of several uniform linear arrays of different elements. For 
example weighted array having element weights [1-2-3-3-2-1] may be 
considered to be made of uniform linear arrays (l-l-l-l-l-l], [l-l-l-l] & 
(1-1] of same element spacing. Some of the weighted elements arrays with 
their components uniform linear arrays are shown in the following Fig.4.7, 
for more illustrations. 


1 . 1 — 2 — 3 — 3 — 2—1 

1 11111 
1111 
1 1 

2 3 - 4 - 5—5 4 - 3 

i 11111 
i 11111 
i 11111 
1111 


3 . 

i 11111 
i 11111 
1111 
i i 

X 1 - 3 - 4 - 3 - 1 

11111 

ill 

ill 


5 . 2 - 4 - 6 - 4 - 2 

11111 

11111 

ill 

ill 


6 Elements + 4 Elements + 
2 Elements 

3x6 Elements + 4 Elements + 
2 Elements 


2x6 Elements + 4 Elements + 
2 Elements 


5 Elements + 2x3 Elements + 
1 Element 


2x5 Elements + 2x3 Elements 
+ 2x1 Elements 


2 - 3 - 4 - 4 - 3 - 2 


Fig. 4.7: Decomposition of weighted array into uniform array 
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Response ofthe array fl-2-3-3-2-11 

To calculate the response R (in dB] of the weighted array (1-2-3-3-2-1] we 
have only to add the amplitude responses ( Ao/ao ) for uniform six elements, 
four elements & two elements and divide by total number of geophones, 
take logarithm and multiply by 20, i.e. 

R = 20 Log 1/12 [{Sin (6 tc d/ X,)/ { Sin(7i d/ X )} 

+{ Sin(4Tt d/X ) / (Sinfic d/ X ]}+{ Sin(27t d/ X )/ { Sin(7t d/ X ]}] 

For reference, the amplitude response of the array (1—2—3—3—2—1] 
have been calculated in the following Table 4.3 and plot array response 
curve has been shown in Fig. 4.8 



Fig. 4.8: Array Response Curve for array 1-2-3-3-2-1 
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Table 4.3 


A/X 

No. of elements 

P 1 +P 2 +P 3 

R=20 Log 
{(Pi+P2+P 3 )/12} 

2 

4 

6 

Ao/ao 
= Pi 

Ao/ao 
= P2 

Ao/ao 
= P 3 

0 . 

2 

4 

6 

12 

0 

0.02 

1.996 

3.961 

5.863 

11.820 

-0.13 

0.04 

1.984 

3.844 

5.462 

11.290 

-0.53 

0.06 

1.965 

3.653 

4.829 

10.447 

-1.20 

0.08 

1.937 

3.395 

4.013 

9.345 

-2.17 

0.10 

1.902 

3.078 

3.078 

8.058 

-3.46 

0.12 

1.860 

2.711 

2.093 

6.664 

-5.12 

0.14 

1.810 

2.307 

1.131 

5.248 

-7.18 

0.16 

1.753 

1.878 

0.260 

3.891 

-9.78 

0.18 

1.689 

1.438 

-0.46 

2.663 

-13.08 

0.20 

1.618 

1.000 

-1.0 

1.618 

-17.40 

0.22 

1.541 

0.58 

-1.32 

0.794 

-23.59 

0.24 

1.458 

0.183 

-1.43 

0.206 

-35.31 

0.26 

1.369 

-0.17 

-1.34 

-0.151 

-38.00 

0.28 

1.275 

-0.47 

-1.09 

-0.299 

-34.39 

0.30 

1.176 

-0.72 

-0.72 

-0.27 

-32.70 

0.32 

1.072 

-0.91 

-0.29 

-0.13 

-38.91 

0.34 

0.964 

-1.03 

-0.14 

0.07 

-44.20 

0.36 

0.852 

-1.08 

0.53 

0.29 

-32.10 

0.38 

0.736 

-1.07 

0.82 

0.49 

-27.74 

0.40 

0.618 

-1.0 

1.0 

0.61 

-25.76 

0.42 

0.497 

-0.87 

1.03 

0.65 

-25.26 

0.44 

0.375 

-0.69 

0.92 

0.59 

-26.03 

0.46 

0.251 

0.48 

0.69 

0.45 

-28.42 

0.48 

0.126 

0.24 

0.36 

0.74 

-24.15 

0.50 

0 

0 

0 

0 

- 


Response of other arrays can also be calculated in similar way as given 
below 

Response of the array 13-4-5-5-4-31 

R = 20 Log 1/24 [{3 Sin (6 n A/X)/ {Sin (ti d /X ]} 

+{ Sin (471 d/ X ]/{Sin [ tzA/X ]}+{ Sin [ 2 tz A/X ]/ { Sin (tc A/X ]}] 
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Response of the array 12-3-4-4-3-21 

R = 20 Log 1/18 [{2 Sin (671 A/X)/ {Sin [n A/ X ]} 

+{ Sin (4rt A/X )/ (Sin [it A/ X ]}+{ Sin (27: A/X )/ { Sin (71 A/X ]}] 
Response of the array fl-3-4-3-11 
R = 20 Log 1/12 [{Sin (571 d^]/ {Sin ( 7 : d/ X ]} 

+ {2 Sin (3 tc A/X)/ {Sin (tc d/ X )}+ {Sin [n A/X)/ {Sin ( 7 : d/ X ]}] 
Response of the array 12-4-6-4-21 
R = 20 Log 1/18 [{2 Sin(5Tt A/X)/ {Sin(7t A/X )} 

+{2Sin(37r A/X)/ {Sin(7c d/ X )}+{2 Sin(7r A/X)/ { Sin(7i d/ X )}] 
Similar expressions can be written for other weighted arrays also. 


4.3.3 RESPONSE OF TAPERED ARRAY OF WEIGHTED SPACING 


(Equal Weight & Unequal Spacing) 

The array in which different elements have different spacing is called 
tapered array of weighted spacing. Response of this type of array can be 
calculated as discussed below. Let us consider four elements uniform linear 
array at uniform spacing'd'. 

0 0 0 0 

< -► 

d 

< -► 

3d 

This array may be considered as combination of two linear arrays each of 
two elements, one with spacing 'd' & other with spacing '3d'. Thus to 
calculate the response of 4 elements uniform array at any point, say A / X = 
0.1, we have simply to add the responses of two elements uniform arrays for 
the point d / X = 0.1 & 0.3. This has been shown in the following Table 4. 3 
and which gives the same result as obtained in the previous calculation for 
response for 4 elements uniform array. 
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Table 4.4 


d//. 

Ao/3o= P1 

For d/X 

(Two element) 

Ao/ao=P2 

For 3d/7 
(Two element) 

P 1 +P 2 

0.10 

1.902 

1.176 

3.078 

0.12 

1.86 

0.852 

2.712 

0.14 

1.81 

0.497 

2.307 

0.16 

1.753 

0.126 

1.879 


In the same way weighted element array 1-2-3-3-2-1 of spacing'd' may be 
considered to be the combination of 6 uniform linear arrays of two elements, 
of these three arrays have spacing'd', two arrays have spacing '3d' and one 
array has spacing '5d' as shown below. 


0 


-* 


0 

0 

-* 


0 0 

0 0 

0 d 0 

◄-► 


3d 

5d 


0 

0 

*• 


0 




Response of the array has been calculated in the following Table 4.5 

Table 4.5 


d//» 

Ao/a 0 =P 1 

For d/X 

Ao/a 0 = P 2 

For 3d/?L 

Ao/ a 0 = P 3 

For 5d/?L 

P1+P2+P3 

FOR TWO ELEMENTS 


3x1.996 

2x1.965 

1x1.902 

11.82 


3x1.984 

2x1.86 

1x1.618 

11.29 

0.06 

3x1.965 

2x1.689 

1x1.176 

10.44 

0.10 

3x1.902 

2x1.176 

0 

8.058 
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These results also agree with the previous results calculated in different way 
(see Table 4.3], Thus we can generalize the method of calculation as 
discussed below. Let us consider an array consisting of N uniform arrays of 
two elements, of which mi number of arrays have equal spacing di, m 2 
number of arrays have equal spacing d 2 , and so on. Then amplitude of the 


combined array will be given by 

Ao = ao [{mi Sin (2 tc di Sin 0 / A ] / (Sin (n d 1 Sin 0 / A ]} 

+{m 2 Sin (2 7i d 2 Sin 0 / A, ] }/{Sin (71 d 2 Sin 0 / A ]} +.] 

Ao = 2ao [{mi Cos (71 di/A)} + {m 2 Cos (tt d 2 /A]} +.-] 

Ao = 2ao [X{ m i Cos(te d i/ A] }] 

If sensitivity of each element is not equal, then 
Ao = 2 [X{ a ;m i Cos(tc di/A)}] 


Since total number of detectors in the array is 2N, normalized amplitude 
response in dB will be given by 

R = 20 Log [Ao/ 2Na 0 ] = 20 Logl/N [Xm j Cos(rc d;/ A]] (4.2ii) 

For an example the response of the following array can be calculated directly 
with the help of the equation 4.2 (ii] 

(Here, N=3, mi=l, m 2 =l, m 3 = 1, di= d 2 /4 = d 3 /8] 

o 2di o 1.5 di o di o 1.5di o 2di o 

4-► 

d : = 4dj 

•4-► 

di = 8di 

4-► 

R= 20 Logl/3 [Cos (71 di/A] + Cos {An di/A] + Cos (871 di/A] ] 


4.4 RESPONSE OF AREAL PATTERN 


The linear arrays are effective for attenuation of those noises, which are 
traveling along the direction of array. In the area of rough topography 
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horizontally traveling noise waves after reflections or scattering may come 
from any direction towards the array. In this case use of areal pattern is 
more effective. One areal arrays is shown in Fig. 4.9. The effective array in a 
given direction can be found by projecting the geophone positions along a 
line in that direction for example, the rectangular array of Fig. 4.9 has effec¬ 
tively a three elements uniform 3-3-3 linear array with the element spacing 
'a'. 



Fig 4.9: Areal Arrays Pattern 

But at 45° to the line, the effective array is 5 elements tapered array 
1-2-3-2-1 with element spacing d = a/V2. 

4.5 RESPONSE OF SOURCE PATTERN 

Theoretically, the 6 elements source arrays produce the same results as 
produced by 6 elements receiver array. A six elements receiver array can be 
replaced by six shot points with same array geometry and shots are fired 
simultaneously with 1/6 th of the charge size distributed in each hole. 
However, we use receiver array more frequently than the source array 
because the cost is less. In difficult areas combination of multiple shots and 
multiple geophones may be used simultaneously. Source arrays are 
generally used when using vibrators and air-guns as seismic source and 
have been discussed in chapter three. 

4.6 ATTENUATION OF RANDOM NOISE 

Output of any detector consists of the signal and noise. Let output g;of the i 
th. detector be, 


gi (t) = Ss (t) + n i (t) 
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Where, Si (t] = Instantaneous amplitude of signal at the position of ith 
detector at time't' 

and m (t]= Instantaneous amplitude of noise at the position of the ith 
detector at time ‘t’ 

The mean square output of the ith detector is the time average, 

[gi 2 ] = (1/T) I gi 2 (t) dt, 

Using above equation. 

[ gi 2 ] =[S, 2 (t) ]+ [n i 2 (t)] +2[S, (t) n , (t)] 

Total output of 'N' detectors in the array at any time will be the sum of the 
instantaneous amplitudes for all detectors. 

gN (t) =£gi (t) 

The mean square output of the array of N detectors is then given by, 

[gs 2 (t)] = Z[S i S j ]+Z[n i n j ]+2Z[S i n, ] 

If the noises & signals are uncorrelated from detectors to detectors in an 
array, 

then [S, nj] = 0 for all i, j and the mean square output of the array is 
[gs 2 (t)] = L [ Si Sj ] + Z [ ni nj ] 

If we make the following assumptions 

• The signal is identical at each detector i.e. S; = S j = S 
For all i & j 

• The noise power is identical at each detector i.e. [ni 2 ] = n 2 
For all i 

• The noise is uncorrelated from detector to detector i.e. [m.nj.] 
= 0 for all i # j 

Then the mean square output of the array is given by, 
g s 2 = N 2 [S 2 ] + N[ n 2 ] 

Mean square signal / Mean square Noise = N 2 [S 2 ] / N [n 2 ] 

Root Mean square signal/ Root Mean square Noise =V N { [S 2 ] 1 / 2 /[ n 2 ]V 2 } 
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Thus the rms signal amplitude to rms noise amplitude ratio of the array is 
improved by factor Vn of the ratio of rms signal amplitude to rms noise 
amplitude of the individual detector. Therefore if large numbers of detectors 
are available, they must be distributed at different places for maximum 
attenuation of random noise. 

4.7 OPTIMIZATION OF ARRAY PARAMETERS 
4.7.1 SMITH’S APPROACH 

Maximum allowable seismometer spacing (Smith—1956 ) is determined by 

d max = v L /f max- (4.3) 

Where Vl = Lowest noise velocity. 

f max = The highest frequency component of reflection 

If A l is the longest noise wavelength which must be kept at least at the 
beginning of the reject region (half power point). 

then, N . d = 0.44 A l 

or N .d = Ai/ 2 (Approx.) 

Where, N = number of elements for uniform array 

d = element spacing 

Since A l should be equal to the shortest wavelength of signal, we have 
Nd — V a / (2 f max ) 

Nd = V / (2 f max Sin 0) 

Nd = [V 2 (t/x)] / (2 f max ) (4.4) 

where, d = Seismometer spacing 

V a = Apparent velocity of signal. 

Sin 0= x/ (v . t) & 0 = Angle of approach for the signal for plane reflector at 
far-offset. 

Smith has shown, in the case of dipping reflector, that 
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Nd = [V 2 (t/x)] / [2 f max {1+ Sin (() ((Vt/X) 2 - Cos2 ()) - Sin 4* )}] ( 4 . 5 ) 

Where, t = Reflection time at far offset 
X = The far Offset distance 
<() = Anticipated dip. 

In general, since the last null point in the response curve occurs at 
wavelength (N-1]/N.d, we must put shortest wavelength noise at or before 
this point, Therefore, we must have 

d < [(N-l] /N ]X sn (4.6) 

where, X sn = Shortest wavelength of noise 

From equations 4.5 & 4.6, suitable valves of d & N may be selected. 

4.7.2 PARR & MAYNE'S APPROACH 

Parr, J.O. & Mayne W.H., Geophysics - 1955, Vol.20, (P539-564) "A new 
method of Pattern Shooting" has suggested a criterion of finding optimized 
weighted array as described below. This array is symmetrical about the 
center of the array and all the secondary peaks in the reject region have 
equal magnitudes. Detailed proof of obtaining optimized array of five 
elements has been discussed below. 

Let the array to be designed be a 2 , ai, 1, ai, a 2 which is symmetrical about the 
center. We have to find the values of weights ai, a 2 and spacing 'do' for the 
array to attenuate waves between shortest wavelength X sn and longest 
wavelength 'X l '. 

a 2 ai 1 ai a 2 

0 0 0 0 0 

A B C D E 

The phase difference for the wave of wavelength X between two successive 
elements will be given by 

0 = 211 d 0 /X (4.7) 

Taking reference in the center, 

Output of element at C = 1 + i x 0 
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Output of element at B = ai(Cos 9 - i sin 9) 

Output of element at A = a 2 (Cos 20 - i sin 20] 

Output of element at D = ai (Cos 0 + i sin 0] 

Output of element at E = a 2 (Cos 20 + i sin 20] 

Vector sum of all elements will be 
1 + 2aiCos 0+ 2 a 2 Cos 20 ; and relative amplitude is 

R = [1 + 2aiCos 0+ 2 a 2 Cos 20 ] / [ 1+2 ai + 2a2 ] ( 4 . 8 ) 

The peak relative amplitude will occur when dR/ dO = o 
i.e aiSin 0 + 2 a 2 Sin. 20 = 0 
Sin 0 (ai+ 4a2Cos9] = 0 

The roots of this equation give the values of 0 for the peak relative amplitude 


and are given below 

Sin 9i =0 i.e. 0i=0orIT (4.9) 

Cos 0i=l When 9i = 0 (4.10) 

Cos 0 2 = -1 When 02 = n (4.11) 

Cos 0 3 = - (ai / 4a 2 ) 

Let us define an arbitrary function 

S = (N-l]a 2 /ai (4.12) 

Where, N = Number of elements 

hence, S = 4a2 / ai for N=5 (4.13) 

and Cos 0 3 = -1 / S (4.14) 

Case-I 


When Cos 9i = l,i.e. 0i = 0 ; From eqn. (4.7) d 0 / ^ = 0 

This peak occurs in the case when energy is received simultaneously at all 
elements. 
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Case-II 

When Cos O 2 = -1 O 2 = n ; From eqn. (4.8] d 0 / X = 0.5 
This peak occurs in the central position of attenuation band. 

Case-Ill 

When Cos O3 = -1/S. there are two values of O3 which can satisfy equation 
[4.14], so there are a pairs of peaks corresponding to this value in the 
attenuation band. 

The broadest attenuation band exists for given relative amplitude when all 
of the peaks in the band are equal in magnitude and alternate in polarity. 
Hence, the value of R as determined from equation [4.8] when Cos 82 =- 1, 
should be equal in magnitude to the value of R when Cos O 3 = -1/S but 
opposite in polarity. Hence, 

[1 + 2ai (-1] + 2 a 2 [+1] ] / [ 1+2 ai + 2 a 2 ] 

= - [1 + 2ai [-1/S] + 2 a 2 (2/ S 2 -l) ] / [ 1+2 ai + 2a 2 ] 

2- 2ai- 2ai/S + 4a2/S 2 =, and since, a2=aiS/4 

2 - 2ai-2ai/S + ai/S = 0 (4.15) 

Thus, al =2 S/ (2S+1) (4.16) 

And a 2 = S 2 /(2(2S+1)} (4.17) 

The value of maximum relative amplitude 'E/ in attenuation 

band can be obtained by substituting Cos O 2 = -1 in equation (4.8) which 
gives 

E r = [1 + 2ai (-1) + 2 a 2 ] / [ 1+2 ai + 2a 2 ] 

Substituting the values of ai & a 2 from equation (4.16) & (4.17), the 
maximum relative amplitude in the band becomes 

E r = (S-l) 2 / ( S 2 + 6S + 1) (4.18) 

To solve, for other values of 0, where the amplitude has the same value as 
these peaks, let us set the value of Er from eqn.(4.18) in place of R in 
equation (4.8) and solve for the values as follows: 



no 


(S-l) 2 /(S 2 +6S+l)=[l+2aiCos0+2a 2 {2Cos 2 0-l}]/[l+2ai+2a 2 ] (4.19) 

Substituting equation (4.16) & (4.17) in the above equation and 
after simplification we get 

S Cos 2 0 +2 Cos 0 -S + 2 = 0 

S Cos 2 0 +S Cos 0 - S Cos 0 +2 Cos 0 - S + 2 = 0 

(Cos 0 + 1) ( S Cos 0-S + 2) = O 

Which gives Cos 0 = -1, as obtained in equation (4.10 ) and 
(S Cos 0 - S + 2)= 0 

Which gives Cos 04 = (S-2) / S (4.20) 

There are two values of 04 between 0 and 2 n that satisfy equation (4.20). 
One value 0 l gives the position of longest wavelength A.l, in the band and 
other value 0s is for the shortest wavelength A.sNin the band. Where, 0s = 2YI- 
0L 

From equation (4.7) 

2II do / 7 -sn = 211- 2fl do/J-L 

d 0 = (^sn A.l) / (^sn + A,l) , this gives the value of element spacing 'do'. 

The value of arbitrary function ‘S’ can be obtained with the help of equation 
(4.20), which gives S = 2/(1-Cos 0 l) Where, 0 l = 211 do /Xl 

By similar method we can design optimized array of any number of 
elements. Results are summarized below 

Three Elements 

0 =1 

ai=S/3 

E r =(S-l)/(S+l) 

Four Elements 

ai=l 

a 2 =S/3 



Ill 


E r =[V{(S-l)/S}] (S-l)/(S+3) 

Five Elements 
A 0 =l 

ai=2S/(2S+l) 
a 2 =S 2 /2 (2S+1] 

E r =(S-l]2/(S 2 +6S+l] 

Six Elements 
ai=l 

a 2 =S/(S+l) 

a 3 =S 2 /5(S+l] 

E r = [V{(S-1)/S}] (S-l] 2 /(S2+10S+5] 

Seven Elements 
C = 3S2+6S+1 
a 0 = 1.0 

ai = 3S.(2S+3]/(2c] 
a 2 3S 2 /c 
a 3 = S 2 /(2c] 

E r = (S-l]/ (6S3+15S2+15S+1] 

Eight Elements 
C = S2+3S+1 
ai= 1.0 
a 2 = S (S+2]/C 
a 3 = S 2 /c 
a 4 = S 3 /7c 

E r = [V{(S-1)/S}] (S-1] 3 /(S 3 +21S 2 + 35 S+7] 
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Nine Elements 
C = 4s 3 +18S 2 +12S+l 
a 0 = 1.0 

ai = 4S (S 2 +4S+2]/C 
a 2 =2S 2 (2S+5]/C 
a 3 = 4S 3 /C 
a 4 = S 4 /2C 

E r =(S-l] 4 /(S 4 +28S 3 +70S 2 +28S+l] 

Ten Elements 
C= S 3 +6S 2 +6S+1 
ai= 1.0 

a 2 = S (3S 2 + 15 S + 10]/3C 
a 3 = S 2 (S + 3]/C 
a 4 = S 3 /C 
a 5 =S 4 / 9c 

E r = [V {(S-l]/S}] (S-1] 4 /(S 4 +26S 3 +126S 2 + 84S+9] 
Eleven Elements 
C = 5S 4 +40S 3 +60S 2 +20S+1 
a 0 = 1.0 

ai=5S (2S 3 + 15S 2 + 20S + 5] /(2C] 

a 2 =5S 2 (S 2 + 6S + 5]/C 

a 3 =5S 3 (2S+ 7]/(2C] 

a 4 =5S 4 /C 

a s =SV(2C) 

E r = (S-l] /(S 3 + 45S 4 + 210 S 3 +210 S 3 + 45 S+l] 
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Example 4.1 

To attenuate the noise having wavelength A,sn= 11.25 m and L=45 m by 6 
elements tapered array with spacing do 

Solution: 

do=^sN W( ^sn+ k)-9m. 

0l= 11x9/45=720 
S= 2/ (11-CosO l ] =2.894 
ai= 1.0 
a 2 = 0.743 

a 3 = 0.342 

E r = 0.0855 (=-21.36 dB] 

Resulting array is 

M -—-► 

9 m 

0 0 0 0 0 
0.342x3 0.743x3 1 x3 1 x3 0.743x3 

Which approximates to the array 1-2-3-3-2-1, when each element weight is 
multiplied by 3 

0 0 0 0 0 0 

1 2 3 3 2 1 

4.7.3 CHEBYSHEV APPROACH 

Mark Holtzman (1963] suggested a most practical and effective optimized 
geophone array called Chebyshev Array for cancellation of surface noise. 
This deals with the arrays in which the reject band peaks are all of equal 
amplitudes. These arrays have the remarkable property that for a given 
absolute upper bound in the reject band, they have the narrowest width of 
pass band lobe as measured between the first nulls. He has claimed that this 


0 

0 342x 3 
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array is more effective than the array suggested by Parr & Mayne and by 
Smith. Procedures suggested by him is as follows 

Calculation of Element Spacing *d' 

Since Asn & Al are the cut off points in the reject band and are symmetrically 
situated about do/A = 0.5 

d 0 / Asn + d 0 / Al = 1 

do = Asn Al/(Asn + Al] 

Calculate Parameter *X' 

X= Sec n do/Xi 

Calculation Of Parameter *m' 

If rejection ratio for Al to Asn be R: 1 i.e. response in the reject band must not 
exceed the value 1/R of the pass band peak. For all practical purposes, R>1 
and hence m>l, then, 

m =Cosh 1 R / Cosh- 1 X 

Calculation Of Number Of Elements 

Numbers of elements in the array 

N =m+l 

Calculation Of Weight Factor 

Weight factors 'a k' is given by 

m/2 

a k = 2/N I Tm [X Cos IT S /N] Tm-2 k [Cos n S /N] 

S=0 

Where, E s =l,ifS = 0 
E s = 2, if S t 0 

and Tm [U] is called Chebyshev Polynomial of first kind defined by 
Tm(U] = Cos(mCos-iU], | U | <1 
= Cos(mCosh iU], |U| >1 
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Which, gives 

To(U) = 1 

Ti(U) = U 

Recursion Formula is 

T m+ i(U) = 2UT m (U]-T m -i(U] 

T 2 (U] = 2U 2 -1 

T 3 (U] = 4U 3 -3U 

T 4 (U] = 8U 4 - 8U 2 +1 

Response Curve 

Complete graph can be plotted from the follow eqation 
Sn (do/A] = | T m [X cos (n do Sin0 /A]] \ 

If 0 = n/ 2 (the angle between direction of wave propagation and normal 
to array direction] 

S N (do/ A) = | T m [X cos (fl d 0 /A]] | 

Example 4.2 

To calculate the values of required number of elements & Array length L for 
getting different degree of attenuation for noises as given below 

(a] AL=50m A SN- 10 m 

(b] AL=50m. AsN=20m 

(c] AL=45m AsN=20m 

Solution: 


A sn =10 m 


Case fAl 
Al =50 m 
do =8.33 m 


X = 1.1547 
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Ratio of 
Amplitudes 

R 

R 

(dB) 

m 

N 

L O) 

10 

20 

5 

6 

41.65 

16 

24 

6 

7 

49.98 

17.78 

25 

7 

8 

58.31 

31.62 

30 

8 

9 

66.64 

56.23 

35 

9 

10 

74.97 

100 

40 

10 

11 

83.33 


Case TbI 

A l = 50 m 
do = 14.3 


A sn -20 m 
X = 1.6057 


Ratio of 
Amplitudes 

R 

R 

(dB) 

m 

N 

L(m) 

10 

20 

3 

4 

42.9 

17.78 

25 

3 

4 

42.9 

31.62 

30 

4 

5 

57.2 

56.23 

35 

5 

6 

71.5 

100 

40 

5 

6 

71.5 


Case fcl 

Al = 45 m Asn = 20 m 
do =13.8 m X = 1.7522 
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Ratio of 
Amplitudes 

R 

R 

(dB] 

m 

N 

L (m) 

10 

20 

3 

4 

41.4 

17.78 

25 

3 

4 

41.4 

25 

28 

3.4 

4 

41.4 

31.62 

30 

4 

5 

55.2 

56.23 

35 

4 

5 

55.2 

100 

40 

5 

6 

69 


Example 4.3 

To calculate the weights of the elements to attenuate the noises having 
wavelengths X l = 45 m, and X sn =20 m. by 25 dB. 

Solution: 

As shown in Ex.4.2 (Case-c], this requires 
do= 13.8 m. 

X=1.7522 

L = 41.4 m N = 4, m=3 

m/2 = 3/2 =1.5 = 1 (for calculation purposes, lowest integer value] 

Weights are given by 

m/2 

a k = 2/N X T m [X Cos Ft S /N] T m - 2k [Cos U S /N] 

S=0 

where, Es = 1, if S = 0 
E s = 2, ifS^O 
Therefore, 

1 

a 0 = %I E s T 3 [ X .CosriS/4]. T 3 [Cos IIS/4] 

S=o 
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=% ITs (X. Cos 0], T 3 (Cos 0] + 2T 3 (X .Cos IT /4].T 3 (Cosn /4)] 
=% [16.2618—5.5026] = 5.3796 = 5(approx.] 
and ai=y 2 [T 3 (X. Cos 0] Ti(Cos 0] +2T 3 (X. CosIT/4] Ti(Cosn/4]] 

=y 2 [16.2618 + 42.8209] 

=29.5414 = 2 30 [approx.] 

Thus after normalization, weight factors ao, ai, ai, ao are 1-6-6-1 

Table - 4.6 


do/A 

Sn 

R(dB] 

do/A 

Sn 

R(dB] 

0 

16.26 

0 

0.28 

2.222 

-17.29 

0.02 

16.145 

-0.06 

0.3 

0.999 

-24.22 

0.04 

15.798 

-0.25 

0.32 

0.494 

-30.35 

0.06 

15.232 

-0.57 

0.34 

-0.126 

-42.21 

0.08 

14.462 

-1.02 

0.36 

-0.577 

-28.99 

0.1 

13.512 

-1.61 

0.38 

-0.862 

-25.16 

0.12 

12.409 

-2.35 

0.4 

-0.989 

-24.31 

0.14 

11.184 

-3.25 

0.42 

-0.976 

-24.43 

0.16 

9.874 

-4.33 

0.44 

-0.843 

-25.7 

0.18 

8.514 

-5.62 

0.46 

-0.616 

-28,42 

0.2 

7.142 

-7.15 

0.48 

-0.325 - 

-33.99 

0.22 

5.793 

-8.96 

0.52 

0.325 

-33.99 

0.24 

4.504 

-11.84 

0.54 

0.616 

-28.42 

0.26 

3.304 

-13.84 

0.56 

0.843 

-25.7 


atA= oo , Response S 4 (d / A] = |T 3 [1.7522]!= 16.26 = + 24.22 dB 





119 


at A = 45 m, Response S 4 (d/ A ] = T 3 [1.7522{ Cos(180xl3.8]} / 45 ] = 1.0 i.e. 
Relative response at = 45 m, with respect to A= go in dB will be given by 

20 Log(l/16.26] =-24.22 dB 

at A = 20 m„ response S 4 (d/ A ] = T 3 [1.7522 {Cos (180x13.8]} /20] = 0.8667 

i.e. Relative response at A = 20 m, with respect to A= go in dB will be given by 
20 Log [0.866/16.26] =-25.47 dB. 

Thus we see that response of the array at Al = 45 m & AsN=20m are 
approximately equal. Complete response curve can be plotted from the 
Table 4.6 and first half of the response curve has been shown in Fig. 4.10. 
The second half of the curve is the mirror image and easily can be 
constructed. 



Fig. 4.10: Response curve for the selected array 












120 


4.7.4 PRACTICAL APPROACH 

In Smith's approach, longest noise wavelength has been kept at 3db 
attenuation point, which corresponds to very small attenuation. In Parr & 
Mayne's and Chebyshev's approach, any amount of noise attenuation can be 
achieved if we don't limit the element spacing and array length. But 
attenuation of highest frequency component of signal by this array has not 
been considered at all. Therefore arrays selected by these criteria may not 
be practically useful. Hence I have considered a practical way of selecting 
array parameter on the consideration of signal and noise both. 

Uniform array Vs weighted array 

As can be seen from the response curve, uniform arrays always have lesser 
effective length than the tapered array for equal amount of attenuation of 
noise. Consequently, due to large effective length tapered array produce 
more attenuation on signal as compared to uniform array. Therefore, as far 
as possible, uniform arrays may be used if noise problem is not very severe. 

Optimum Value of (X/Z) 

First of all maximum allowable value of (X/Z] is calculated on the basis of 
shortest wavelength component of signal required to fulfill the geological 
objectives and prominent noise wavelength which is interfering with the 
signal at the level of main zone of interest. Attenuation of desired shortest 
wavelength signal must not be more than 3 dB (half power point] and that of 
prominent noise must be more than or equal to 20 dB. As shown in Fig. 2.5, 
3 dB attenuation occurs at L / A sx = 0.44 and 20 dB attenuation at L / A n= 
0.91 for any uniform array of elements more than 3. Where, A sx is the 
apparent shortest wavelength of signal and An is the wavelength of 
prominent noise and 'L' is effective array length. Thus we have from eqn. 2.6 
and 2.7as follows. 


L / A sx < 0.44 

L/A n > 0.91 

And for limiting values 

A N /A SX =0.44/0.91= 0.48 


(4.21) 

(4.22) 
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A n /A s = 0.48 Cosec 0, since A sx = A s Cosec 0 

Substituting the value of Cosec 0 (for dipping reflector, Ch.2), we get 
An/ 0.48A s = 1+ [4Z2.Cos20/ (X±2Z Sin0)2]% 

After simplification this becomes 


X/Z = 2 Cos 0/[(2.O8 An/As) 2 -1]V4 ± 2Sin 0 


(4.23) 


Here, (+) for up dip shooting and (-) for down dip shooting. 

Now, value of (X/Z) is calculated, on the basis of the value of A n,As and 0 (dip 
angle). If (X/Z) calculated, comes out to be equal to or greater than the value 
previously decided in Chapter-2, then there is no problem, we can keep the 
previous value of X/Z. If it is less than the previous value, then we must keep 
less far offset as calculated here, to attenuate the noise and preserve the 
signal at required level. Otherwise we have to satisfy with, either lesser 
attenuation of noise or greater attenuation of signal. So, a compromise for 
(X/Z) value is done at this level, if required. 

Effective array length 'L' 

After deciding (X/Z) value, the maximum effective array length is calculated 
on the basis of following relations: 

Putting, Asx = 0.44 As Cosec 0 in eqn. 4.21 

L = 0.44 Asx = 0.44 As Cosec 0 

Substituting the value of Cosec 0 for different cases, we get the follow result 

Up-dip 


L = 0.44 As [1+ (4Z2Cos2 0/ (X-2Z Sin 0) 2 )] % 


(4.24) 


Down-dip 


L = 0.44 As [1+ (4Z 2 Cos 20/ (X+2Z Sin 0) 2 }] % 

NO-DIP 


(4.25) 


L = 0.44 A s [1+ (4Z 2 / (X) 2 )] % 


(4.26) 
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Number of elements 'N' and element spacing' d' 

If value of'L' calculated from the above relation is more than the prominent 
noise wavelength A n, then we may keep for putting X n at the first notch 

A N =Nd for A n <L (4.27) 

Otherwise, L = Nd for A n > L (4.28) 

to keep A Nat 20 dB attenuation point. 

N & d may be varied, such that N > 4 and d < geophone spacing in the string 
and base-length closest to group-interval, if equal values of base-length and 
group interval are not possible. 

No. of elements, N= L/(L-Lb) (4.29) 

Where, Lb = desired base length of the array, preferably equal to group 
interval. 

If equation (4.29) does not provide suitable value of N, then any other 
suitable combination of N & d may be tried satisfying the equations (4.27) or 
(4.28). and their effectiveness can be tested by fold-back experiment and (F- 
K) diagram. 

Example 4.4 

In an area, An = 36 m, As = 40 m, Z = 2500 m, 

Design the optimum array. (Assuming 0 = 0) 

Solution: 

As calculated in example 2.8, X = 3000 m, AX = 30 m, 

(i) This gives, X/Z = 1.26, according to eqn. (4.23), i.e. X=3150m, which 
comes out more than the previous value, X = 3000 m (as given). 
Therefore we may keep, X= 3000m safely. 

(ii) Now, X = 3000 m, Z = 2500 m. 

Therefore, L = 34 m, according to eqn.4.26 

(iii) Since, A n > L, 

We put, L = Nd (eqn.4.28) 
i.e. Nd = 34 
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(iv) If we put base-length = group interval = 30m (as given] 
No. of elements, N = L /(L-Lb) = 34/4 = 8.5 = 8 (Say). 

(v) Element spacing'd' = L /N= 34/ 8 = 4.25 m. 


4.8 FOLD-BACK EXPERIMENT AND ANALYSIS OF THE RESULT 

As discussed above, two possible optimum arrays are decided, based on the 
signal and noise characteristics and their effectiveness is tested with respect 
to the bunched geophones in real field conditions. Layout of the array is 
shown below. If'n ‘ is the numbers of active channels, channels 1 to (n/3) 
are in the first row and geophones are planted with parameters of first array 
pattern. In the second row at a distance of a few m, are channels (n/3+1) to 
(2n/3) and geophones are planted with bunching of geophones, and in the 
third row are channels. (2n/3+l) to (n) with parameters of second array 
pattern. Group interval AX in each row is the same and equal to the group 
interval to be kept in production work. First shot with optimum parameters 
is taken at an offset, equal to near offset to be used in production work. 
Second shot at an offset equal to (nAX/3) from the first shot point, and third 
shot again at an offset of (n AX/3) from the second shot point are taken. For 
96 channels recording, 30m group interval and initial near offset 150 m, first 
record shows the behavior of arrays from 150 m to 1080 m. Second record 
shows the behavior of arrays from 1110 to 2040m. Third record shows the 
behavior of arrays from 2070m to 3000m. 3000 m is the far offset distance 
to be kept in production work. 

Channels having the same offset in the three groups of detectors are 
compared for each offset. S/N and peak frequency of the signal is calculated 
at the zone of interest, from the refloat monitor record. Amplitude spectrum 
and auto-correlation function can also be obtained at the zone of interest at 
several offsets with the help of computer. In the field lay-out as shown 
below, channels 1 to 32 are having array 3-4- 5-5-4-3 with geophone spacing 
6m; channels 33 to 64 are having bunched geophones and channels 65 to 96 
are having array 4-4-4-4-4-4 with geophone spacing 6 m. Fold-Back layout is 
shown below in Fig. 4.11 and field monitor record in Fig. 4.12. 
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Fig. 4.12: Fold Back Field Record 

Frequency analysis has been done at four offsets and plot of frequency 
spectrum are shown in Fig.4.13 (offset 1080 m] and Fig.4.14 (offset 150 m]. 
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Channels 1, 64 & 65 correspond to an offset = 1080 m, channels 16, 49 & 80 
correspond to an offset = 600 m, channels 32, 33 & 96 corresponds to an 
offset = 150m and channels 13, 52 and 77 correspond to an offset = 510 m. 

Table 4.7 



Offset 1080 m. 

Offset 600 m. 

Offset 150 m 

Offset 510 m. 


Channe 

S 

Channels 

Channels 

C 

hannels 

FREQUENCY 

1 

64 

65 

16 

49 

80 

32 

33 

96 

13 

57 

77 

NORMALISED 

NORMALISED 

NORMALISED 

NORMALISED 


AMPLITUDE 

AMPLITUDE 

AMPLITUDE 

AMPLITUDE 

0 

4 

3 

1 

1 

1 

0 

2 

3 

2 

8 

4 

10 

2 

4 

4 

1 

2 

1 

2 

3 

3 

2 

8 

4 

10 

4 

4 

4 

1 

4 

2 

4 

4 

5 

1 

9 

4 

11 

6 

4 

4 

2 

6 

3 

7 

6 

6 

1 

9 

4 

11 

8 

5 

4 

3 

8 

6 

8 

8 

12 

3 

13 

7 

20 

10 

6 

4 

2 

8 

9 

9 

12 

0.4 

5 

18 

8 

11 

12 

11 

10 

15 

39 

22 

41 

38 

44 

34 

59 

34 

62 

14 

43 

48 

34 

100 

77 

48 

31 

94 

38 

100 

100 

92 

16 

61 

45 

43 

78 

100 

71 

80 

62 

60 

40 

49 

64 

18 

8 

40 

15 

92 

78 

100 

48 

54 

52 

89 

51 

91 

20 

22 

19 

34 

48 

27 

45 

100 

100 

86 

40 

20 

60 

22 

100 

100 

100 

42 

17 

26 

41 

37 

40 

39 

25 

45 

24 

12 

12 

18 

21 

10 

20 

14 

15 

31 

38 

16 

27 

26 

41 

27 

47 

74 

32 

50 

81 

66 

100 

74 

31 

100 

28 

33 

22 

43 

90 

40 

67 

69 

55 

79 

56 

26 

84 

30 

16 

5 

17 

14 

4 

11 

37 

19 

30 

11 

8 

2 

32 

10 

5 

12 

18 

6 

17 

15 

19 

26 

31 

5 

47 

34 

16 

5 

20 

7 

3 

7 

21 

19 

26 

20 

14 

42 

36 

17 

8 

25 

30 

14 

19 

29 

24 

34 

30 

10 

29 

38 

13 

11 

18 

18 

6 

12 

30 

15 

29 

24 

7 

31 

40 

2 

2 

6 

14 

6 

10 

8 

10 

15 

8 

4 

7 

42 

1 

2 

2 

20 

5 

15 

10 

8 

11 

22 

10 

28 

44 

4 

3 

4 

20 

9 

18 

5 

2 

6 

21 

6 

26 

46 

1 

1 

1 

15 

7 

15 

11 

9 

14 

6 

1 

8 

48 

3 

3 

4 

10 

6 

10 

5 

6 

9 

8 

2 

6 

50 

3 

2 

3 

7 

4 

8 

2 

1 

4 

5 

0 

4 


Table 4.7 shows the Amplitude Spectra of Fold Back Record. Channels 13, 
16 and 57 are having peak frequency 14 Hz (Table 4.7], which is the 
frequency of prominent noise. At all the four offsets, channels 65 to 96 
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(array 4-4-4-4-4-4] are showing higher peak frequencies as compared to 
others. 


- Array 3-4-5-5-4-3 

- Bunched 

- Array 4-44-4-44 



Frequency 


Fig 4.13 Amplitude spectrum of trace at an offset of 1080 m 


Channels 1 to 32 : Array 3,4,5 - 5, 4,3. 

Channels 33 to 64 : Bunched geophones. 

Channels 65 to 96: Array 4,4,4 - 4,4,4 
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4.9 (F-K) DIAGRAM 

(f-k ] diagram is very useful tool for design and performance evaluation of 
an array. This provides a direct indication of how well noises are attenuated 
and what affects the array has on desired reflection events. To draw (f- k] 
diagram, each (f, k] pair both for noise and signal is plotted on (f-k] plane as 
a point ( k is apparent wave number, f is frequency] and these points are 
joined with the origin by a straight line. Slope of each straight line gives the 
corresponding apparent velocity. A practical example of (f-k] diagram is 
shown in Fig. 4.15 and array response curve in Fig.4.16. Signal & noise have 
following characteristics: 
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Wave Number(k) 


Fig.4.15: f-k diagram 



Fig.4.16: Array response curve for 6-element uniform array with 
6m spacing (A is signal and a, b & c are noises] 

Characteristics of Signal 

Let depth of target horizon be Z =2500m, average velocity to the target 
horizon = 2400m/sec and far offset X=3000m. If as per objective, resolution 
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of 10 m thick bed is required, then shortest wavelength component of signal 
required to be recorded is 40 m and corresponding frequency would be fi 
=V/A = 60Hz. The required apparent wavelength of the signal would be 

A x = A /sin 0 = 78 m : where, 0 = tan^X/ZZ] = 31° Therefore, Ki = 0.013 
cycle/m. 

This (f, k] pair corresponds to point (A] in the [f-k] diagram. 

Characteristics of Noise 

According to figures 4.3 & 4.4, the three groups of noises have the following 
characteristics 


(a] 

f 2 = 14Hz, 

V 2 = 440 m/Sec., A 2 = 31 m, 

K 2 = 0.032 

Cb) 

f 3 = 9 Hz, 

V 3 = 370 m/Sec., A 3 = 41 m., 

K 3 = 0.024 

(c) 

U = 5 Hz, 

V 4 = 295 m/sec., A 4 = 59 m., 

K 4 = 0.017 

These are shown 

as points [a], (b] and [c] in the [f-k] diagram. 


Below this diagram, response curve of one suitable array 4 - 4 - 4 - 4 - 4-4 with 6 
m. spacing is drawn. Attenuation of noise (a] is around 19 dB and that of 
noise [b] is about 17.5 dB by this array. Noise (c] is about 7.5 dB and array is 
not effective for this noise. However, as the frequency is very low, noise c 
may be suppressed with the help of low cut filter of the recording device. 
The required signal, which corresponds to point (A] in the (f-k] diagram, is 
well within the Pass-band of the array response curve. Therefore the array 
4-4-4-4-4-4, for the above mentioned noise and signal is a suitable one. 


4.10 OTHER CONSIDERATIONS 

Deviations from the optimum parameters in respect of following may 
seriously affect the attenuation characteristics of the array. 

[a] If geophone plantations are not uniform i.e. either the coupling with the 
ground is different or geophones are inclined at different angles with the 
vertical, then even if the numbers of geophones are equal at each element 
positions, their sensitivity will not be equal. And, therefore, uniform array 
will change into tapered, and position of notches & cut off point will be 
shifted. 
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(b] If geophones are not planted exactly at the same distance as planned, 
even if plantations are uniform, there will be tapering in spacing. In this case 
also characteristics of the array will be changed considerably. 

(c] If there is elevation changes within the group, there will be phase shift of 
the reflected wave coming nearly vertical to the ground within the array 
between element to element, and resulting wavelet will have less amplitude 
and frequency. 

Thus proper geophone plantation is very essential to get desired result from 
geophone array. 
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CHAPTER - 5 


RECORDING SYSTEM AND QUALITY CONTROL 


5.0 INTRODUCTION 

O nce the parameters for data acquisition are decided on the basis of 
given objective of the survey, the prime goal of the survey team should 
be to adhere to those parameters. But in real field condition some times 
slight deviations from the set parameter may occur due to logistics or 
economic reasons. Therefore it is always necessary to put a limits on the 
amount of deviations, so that fulfillment of objective may not be affected 
much. 


5.1 DAILY TEST OF RECORDING SYSTEM 

Several daily tests are performed and results are analyzed. There are pre¬ 
defined limits for these tests, which should not be crossed since otherwise 
the required quality cannot be satisfied. For faithful recording of the seismic 
signal, the amplitude-frequency response of the complete instrument system 
should be flat over the frequency range of the signal and phase-frequency 
response must be a straight line passing through the origin or a multiple of 
2fl on phase axis. But these two requirements are very difficult to achieve in 
practice. Apart from these, all the channels in the instrument should have 
the same gain, the same amplitude - frequency and phase-frequency 
response and the same time alignment of the channels. Behavior of the all 
parts of the instrument should be linear (i.e. independent of magnitude & 
polarity of the signal]. But the linearity of the system is obeyed only up to a 
certain maximum amplitude of the input signal. When amplitude of input 
signal exceeds this limit, some new frequencies (integral multiples of the 
frequency of input signal] are produced in the output and this effect is called 
harmonic distortion. 
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Harmonic Distortion Test 

Harmonic distortion is produced by lack of linearity in the system. It is 
measured by applying a digitally synthesized sine wave with a distortion 
less than 0.01 % to the input. The out is then analyzed to measure the 
amplitude of the harmonics. There are two ways of measuring non-linearity 
in terms of harmonic distortion in percentage i.e. 

Percentage first harmonic distortion = (A2/A0) x 100% 

Where Ao is the amplitude of input signal of frequency ‘f ' and A2 is the 
amplitude of output signal of frequency ‘2f'(the first harmonics]. 

The second is the total harmonic distortion i.e. 

Total harmonic distortion =(A2 2 + A3 2 + A4 2 +. 3 1 / 2 / A 0 x 100% 

Where, A3, A4, . are the amplitudes of 2 nd , 3 rd ,.harmonics and so on. 

If there are two or more frequencies in the input signal then inter¬ 
modulation distortion is developed in the output. If we consider just two 
frequencies fi & f 2 then, the frequencies having the sum and difference of the 
frequency fi and f 2 and all their harmonics i.e. fi + f 2 , fi - f 2 , 2fi + f 2 , fi+ 2 f 2 , 2fi - 
f 2 , fi - 2f 2 etc. will appear in the output. Total harmonic distortion should not 
be higher than 0.025 % for amplitudes between 0 and -40 dB. Harmonic 
distortion varies with frequency and generally more severe at low 
frequencies. Test should therefore be executed with at least three different 
input signals; two signals with a frequency near the high and low cut filter 
edge and one signal with a frequency in the mid-band. 

In the process of deconvolution, to improve geological resolution by 
changing earth filtered seismic pulse back into a simple spike, the amplitude 
spectrum of the seismic signal becomes broadened i.e. low as well as high 
frequencies relative to the dominant frequency is boosted in the spectrum. 
Then if the distortion product is present in the signal, this will also be 
increased and will produce undesirable effect on the signal. Sometimes 
amplitudes of surface waves may be many times greater than the desired 
reflection signal, but their frequencies are below the dominant frequency of 
the reflection. If we allow the harmonic distortion to occur without using 
filtering, the distortion products from the surface waves may be greater than 
the desired reflection signal within the signal spectrum. 





133 


Amplifier Impulse Response test 

For this test input is the Fourier transform of an impulse so that 
performance of the system can be viewed at all frequencies. The output of a 
Fourier transformed 1 ms input pulse to the amplifier is measured for each 
channel. By inspection of the output amplitudes in the frequency domain the 
amplifier's performance can be verified. 

Amplifier rms noise and dc offset test 

Noise is always present in electrical system but should be as low as possible 
since it dictates the minimum observable output. The strongest component 
generating the instrument noise is the thermal noise, resulting from the 
agitation of the electrons in the conductors under the influence of 
temperature. The random motion of these electrons is the equivalent of 
small currents, which develop a voltage at the terminals of the conductor. If 
input of the system is Ym,the gain G and the experiment is conducted ‘ r ' 
times, we obtain the I th output as 

Yw=G (Y in + Ni] 

Where Ni is the value of the system noise for the i th. experiment. Since the 
average value of the of Ni is zero, the estimated value of the output is 

Y 0 ut = G Ym, therefore 
Y’out - Yout = GNi 
(Yi 0 ut-Yout) 2 = (GNj) 2 


N; 2 = (YW- Yout) 2 / G 2 


RMS noise =[ £( Ni 2 / r)]i/2 = [£ (Y out - Y out ) 2 / r]V 2 / G 

If the input signal to a system is zero, one would expect an output signal of 
zero amplitude. If the measured amplitude is not zero it is called the DC 
offset. If DC offset is +Vo, then the amplitude of the largest acceptable signal 
becomes V m -Vo, where V m denotes the largest acceptable voltage in the 
absence of offset. 
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Channel Gain Accuracy Test 

Before digitization of signal, fixed gain is applied to all the channels. This 
gain should be equal in all channels so that lateral variation in amplitude 
along an event may be truly recorded. This test is conducted by feeding the 
same signal to all inputs. The channel-to-channel deviation errors should not 
be higher than 0.5 %. 


5.2 IMPORTANCE OF LOW-CUT FILTER 

Geophone introduces phase distortion of the wave below the natural 
frequency, so signal at these frequency is distorted (for a 10 Hz geophone 
with 0.7 critical damping the lag varies by only 4 ms. between 10 Hz and 100 
Hz but 45 ms. between 4 Hz and 15 Hz. Therefore, minimum value of Low 
Cut should be equal to geophone natural frequency. If high amplitude low 
frequency ground roll be a serious problem and cannot be removed by other 
methods (by processing or by the use of array], a higher low cut filter may 
be desirable to reject such noise. But, generally there is choice between 
reducing ground roll by using long arrays (which may attenuate high 
frequency component of signal] and a higher low cut filter (which may 
attenuate low frequency component of signal]. The choice will depend on 
the target depth and required resolution. For the best resolution, the higher 
low-cut filter is most desirable. Seismic noise, in general are of low 
frequency and their amplitudes are large as compared to amplitudes of the 
signal of high frequencies because seismic sources have more low frequency 
energy than high frequency energy and this spectral imbalance is further 
accentuated by earth, attenuating higher frequency more severely. The 
higher frequencies received at the geophone can be of such small amplitude 
compared to the low frequencies that they do not register in AID conversion. 
They are lost in the system noise, or they are smaller than the voltage value 
represented by the least significant bit of the A/D converter. Application of 
low cut filter can improve this situation by attenuating the low frequencies 
to a level where their magnitude is comparable to the magnitude of the high 
frequencies. Since earth behaves like high cut filter, let us suppose that the 
slope of earth filter is 24 dB/Octave. If the instrument records 14 bits of data 
and peak frequency recorded is 20 Hz. The gain of the amplifier will be set 
almost to saturation (i.e. all 14 bits occupied] at 20 Hz. At 40 Hz the signal 
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will be attenuated by 24 dB with respect to 20 Hz frequency (i.e. loss of 4 
bits] by the earth and the instrument will record only 10 significant bits of 
that frequency. At 80 Hz the signal will be attenuated by 48 dB (i.e. loss of 8 
bits] and only 6 bits of data will be recorded of that frequency component 
At 160 Hz the signal will be attenuated by 72 dB (i.e. loss 12 bits] and only 2 
bits of data will be recorded of that frequency component, which is very 
insignificant. Thus in this situation it is impossible to record frequencies 
greater than 80 Hz at significant levels with low-cut out. For recording of 
frequencies greater than 80 Hz at significant level we must apply low-cut 
filter before AID converter to counter the action of earth high cut filtering 
action. For instance if, low cut filter with cut-off frequency 80 Hz and slope 
24 dB/Octave be applied before A/D conversion the amplitude of 20 Hz 
signal will be 48 dB down, amplitude of 40 Hz signal will also be 48 dB down 
(24 db by instrument filter + 24 dB by earth], amplitude of 80 Hz signal will 
also be 48 dB (by earth filter only]. Therefore, the amplitude of 20, 40 6 80 
Hz frequency components will all be the same and instrument gain now 
saturate at his level. The signal will be attenuated by 24 dB at 160 Hz and 48 
dB at 320 Hz and thus recording of frequencies higher than 100 Hz will be 
possible (assuming our source generated them]. Steeples and Knapp, 
(1982], made an experiment. They first recorded data with low cut out 
(broad-band-recording] and latter applying low cut filter in the identical 
field conditions. The data recorded with broadband recordings were found 
narrower band than the data recorded with the application of low cut filter 
as shown in Fig. 5.1. In practice it is desirable to have a signal with relatively 
flat amplitude spectrum over two to three octaves of frequencies. Otherwise 
the data may show "ringingness" on the record section resulting in poor 
resolution. A frequency band of 10 to 40 Hz and 40 to 160 Hz, have the same 
band-width of two octave but the pulse of band 40 to 160 Hz have greater 
resolution because it has narrow pulse width. Thus for high- resolution work 
larger band-width and higher frequency content is desirable. 

As discussed earlier system becomes non-linear in behavior and produce 
new frequencies in the output if the input is larger than the specified. These 
new frequencies are integral multiple of the input frequencies. Some of this 
may lie in the pass band and this cannot be removed by applying low-cut 
filter at the time of processing. This effect must be removed before 
digitization of the signal in the field itself. Thus if the signal to noise ratio of 
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seismic event is comparable to the recorder signal to noise ratio, noise must 
be attenuated before recording, to ensure that the desired reflection signal is 
recorded safely above the distortion products of the interfering noises and 
safely above the system noise. If entire seismic system were linear, this 
could as well be done in the final processing as in the field. An A/D converter 
is designed to accept a certain maximum voltage depending upon the digital 
word length. The maximum detectable voltage is the level at which the A/D 
converter becomes saturated and its output exceeds the largest possible 
digital word length. Signals above the maximum level are truncated. Clipping 
circuits limit voltage input to the A/D converter. Saturated or clipped signals 
appear as square rather than sinusoids. 



Fig.5.1: Normalized amplitude Spectra with and without filter 
(Courtesy - Steeples and Knapp, 1982] 



137 


5.3 HIGH-CUT FILTER AND SAMPLE RATE 

Sample rate is related with high cut filter. The high cut filter should be 
enough that the highest frequency required in the target reflection i.e. the 
highest frequency needed for getting vertical resolution, or the maximum 
frequency of signal which can be recorded at a given depth, whichever is 
less, is not attenuated. Maximum frequency of the signal vary from area to 
area, but an empirical guide line (Denhan, 1981] is 

f max = 150/t 

Where't' is two way reflection time in seconds and f max is the maximum 
useable frequency. Once the choice of high cut (or anti, alias filter) is fixed, it 
should be compatible with the system sample rate. A minimum of two 
samples per cycle is necessary for avoiding aliasing. If sampling interval be t, 
then frequency '1/ t’ is called sampling frequency and half of this frequency 
is Nyquist frequency (fN). If signal is having frequency greater than the 
Nyquist frequency say this is (fN + f) then the frequency of the recorded data 
will appear to have lower frequency equal to (fN - f). This low frequency 
effect in the output, when the input is having frequency greater than the 
Nyquist, is called aliasing. As a rule of thumb sampling interval is taken as 

t = % fmax 


5.4 DEFINITION OF BAD RECORDS/CHANNELS 

Recorded data are considered to be bad if 

(i) Records/Channels are having 

(a) No seismic data recorded or incomplete record 

(b) No navigation data recorded or incomplete record 

(c) Loss of time break 

(d) Tape parity error 

(e) Channel is dead 

(f) Channel is spiky, leaky or wild 

(g) Sensitivity deviation more than 6 dB from surrounding channels 

(h) Reverse polarity 

(i) Number of streamer compasses do not meet 
performance specifications (marine data) 
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(j) Streamer depth out of specification (marine data] 

[ii] Source does not meet performance specification 

[a] No shot fired 

[b] Source timing or synchronization not proper 

[c] Source separations not proper 

[d] Source volume not proper 

[e] Air pressure not proper 

[f] Source depth not proper 

[g] Failure of source depth indicator 

[h] Malfunctioning of source controller 

[i] Auto fire of gun 

[iii] Recording device does not meet performance specification 

[a] Recording parameters not proper 

[b] Recording instrument test results not proper 

[c] QC system not operating 

[iv] Streamer does not meet performance specification (marine] 

(a] Streamer feathering exceeds performance specifications 

(b] Streamer noise levels exceeds performance specifications 

(c] Number ofbad channels do not meet performance specifications 

[v] Positioning data does not meet performance specification (marine] 

(a] Vessel primary positioning data out of specification specifications 

(b] No gyro information 

(c] Fathometer not operating 

(d] Number of streamer compasses do not meet performance 
specifications 

(e] Streamer depth out of specification 


5.5 3-D MARINE SEISMIC ACQUISITION QC MEASURES 

3-D quality (QC] system should have real-time monitoring of sub-surface 
midpoint coverage for 3D seismic surveys. The hardware and software shall 
provide continuous monitoring of vessel positioning and a continuous 
recording of all auxiliary data, which shall include but not limited to, items 
such as Doppler- sonar and cable compasses. 3D QC system should 
continuously receive required data for requisite geometrical sub-surface and 
surface position and midpoint calculations and data logging. 
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The 3D QC system shall perform calculations of subsurface coverage on an 
orthogonal grid system using the Universal Transverse Mercator (UTM] 
projection. After definition of the spheroid parameters and the projection of 
the grid, the systems shall monitor the midpoint multiplicity in each bin for 
up to four user defined cable segments. Complete grid coverage shall be kept 
on disc with a backup on the magnetic tape or other media. Quality control 
displays shall be monitored in real time and in order to maximize the data 
acquisition efficiency shall include the following on line real-time displays or 
printouts. 

a] Midpoint coverage displayed on a cathode ray tube (CRT] monitor 
device. 

b] Total and partial grid display on plotter or CRT. 

c] Compass statistics printouts, which shall be available at the end of 
each line. These statistics shall include the individual mean and 
standard deviations. In addition, for each line the average cable 
angle and the midpoint angle for entire line shall be printed in. the 
statistical list. 

d] Automatic indication of bad compass values whilst shooting. 

e] Number of midpoint per cable segment on the pre-plotted shooting 
line. 

f] Coverage prediction utility for planning subsequent lines or infill 
lines. 

g] CRT displays shall be output to a thermal printer for hardcopy as 
required. 

h] Edit facility should be available to remove shot-points that have poor 
quality data e.g. bad navigation, misfire and files not recorded. 

i] The system shall allow for dividing the cable into a maximum of 4 
segments and furthermore shall have the capability of verifying that 
the offset distribution specifications have been achieved. For 
example, the required presence of full fold coverage and 
corresponding allowable minimum number of CDP traces falling in 
each bin can be as follows : 
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Segment 

Offset 

(mts) 

Ideal No of 

traces 

Minimum No 

of traces 

Minimum 

Full fold 

coverage 

1 

0 -600 

12 

10 

80 

2 

600-1200 

12 

10 

80 

3 

1200 -1800 

12 

9 

75 

4 

1800-2400 

12 

7 

60 


j] Each bin should maintain full set of trace offsets. Same offset traces 
in each CDP bin will not be considered or fold purpose, If any 
segment does not meet the coverage specified above for more than 
25 consecutive pops or 15% of total line re-shoot / infill has to be 
done. 

k] The system shall be capable of providing a practical display in real¬ 
time of the midpoint distribution of reflection points from each cable 
segment. This shall be referred to when shooting infill lines. 

l] All midpoint coverage data shall be recorded on magnetic tape in 
case of a failure of the data base storage system. 

m] The gyro alignment shall be checked by center lining before the start 
of the survey. 

Based on above considerations, standard quality norms have been fixed, 
which are adhered by the industries. 
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